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^ ABSTRACT 

^ We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extra- 

galactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS 
Evolutionary Probe (PEP) and GOODS-//eri'c/?e/ key programmes. We describe data reduction and the construction of images and 
catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3cr depths of 0.9, 0.6 and 1.3 mjy at 70, 100 and 160 yum, 
respectively, and resolve ~ 75% of the cosmic infrared background at 100 fim and 160 /im into individually detected sources. We use 
these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths, and to determine 
the infrared luminosity function of galaxies down to Lir = 10" Lq at z ~ 1 and Ljr = 10'^ Lq at z ~ 2, respectively. For the infrared 
luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate 
infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (> 3cr) are 
now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide 
a powerful new tool for studying galaxy evolution over a broad range of redshifts. 
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1. introduction 

The detection of a cosmic infrared background, as energetic as 
the optical/near-infrared background (' Puget et al.|1996 Hauser 



|et al.|199 8 ), has revealed the importance of the energy absorbed 
by the dust in galaxies and re-emitted at mid- to far-infrared 
wavelengths. From this finding it became clear that a complete 
census on the formation and evolution of galaxies could not be 
obtained without accounting for this dust emission. Since then, 
many studies have confirmed the importance of dust emission us- 
ing individual detections of galaxies from infrared facilities such 
as the Infrared Space Observatory, ISO, or the Spitzer Space 
Telescope. However, because of their relatively small onboard 
optics, the far-infrared capabilities of these observatories were 
strongly limited by source confusion. At the time, far-infrared 
studies were restricted to the analysis of local galaxies or to the 
analysis of rare high-redshift very luminous galaxies. As a con- 
sequence, only a small fraction of the cosmic far-infrared (i.e., 
A > 40 /vm) background was resolved into individual objects, 
so our knowledge of the high-redshift Universe at far-infrared 
wavelengths was very incomplete. 

Tha nks to the advent of t he PACS (Poglitsch etal.|2010[ ) and 
SPIRE ( [GrifRn et al.||2010|l instruments onboard the Herschel 
Space Observatory (Pilbratt et al. 2010), this limitation has been 
largely overcome. Indeed, using the relatively high spatial reso- 
lution (provided by a 3.5 m mirror) and sensitivity of Herschel, 
deep extragalactic surveys can be pursued, thereby resolving a 
large fraction of the cosmic far-infrared (i.e., ~58% and ~74% 
with PACS at 100 and 160 //m, respectively; |Berta et "aL]|20I0 



2011 1 and submillimetre (i.e., -15% with SPIRE at 250 yum; 



Oliver et al.|2010| l background into individually detected galax- 



ies. From these observations, one can study the origin and nature 
of the cosmic infrared background through, e.g., the determina- 
tion of the infrared luminosity functions of galaxies (e.g., |Grup^ 



of their spectral energy distributions (e.g., Hwang et al. 2010; El- 


baz et al. 2010, 2011; Nordon et al. 


2010 2012'; Magnellietal.| 


2010,,,2012,.Symeonidis et al.|2013 


Berta et al. ,2013(1. All these 



studies point towards the diversity and redshift evolution, both in 
term of numbers and properties, of the infrared luminous galaxy 
population. Relatively rare in the local Universe, the infrared 
luminous galaxies dominate the cosmic star-formation history 
at z > 1 (e.g., |Gruppioni et al.||2010 20131, and their physical 
properties differ significantly from those of their local counter- 
parts (e.g., lElbaz et al.|[20n] |Wuyts et al.|[20n] [Nordon et al^ 
[2012) 1. 

In this context, we study here the infrared luminous galaxy 
population further by deriving the PACS numbers counts and in- 
frared luminosity functions down to unprecedented depths using 
the combination of the two main extragalactic surveys designed 
to take advantage of the f ull PACS cap abilities: the PACS Evolu- 
tionary Probe (PEfQ L utz et al.|201l| guaranteed ti me key pro- 
gramm e; and the GOODS-Herschel (GOODS-FQ |Elbaz et al. 
[20111) open time key programme. The PEP survey is structured 



as a "wedding cake" (i.e., with large area shallow images and 
smaller deep images) and includes many widely studied blank 
and lensed extragalactic fields, such as the Great Observatories 
Origins Deep Surveys North (GOODS-N) and South (GOODS- 

* Based on observations carried out by the Herschel space observa- 
tory. Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

' http : //www . mpe . mpg . de/ ir/Research/PEP 

^ http : //hedam. oamp . fr/GOODS-Herschel 



S) fields and the cosmological evolution survey (COSMOS). The 
GOODS-H survey only focuses on the GOODS fields, but using 
very deep observations, close to the Herschel confusion limit. 
The extensive observations of the GOODS fields made by the 
PEP and GOODS-H surveys is explained by the availability of 
a deep multi-wavelength database incl uding X-ray (A lexander] 
[eraL[[2003| [Xue et al.[[20TT]i, optical ([Giavahsco et al.|i2004l i7 
near- nfrared (CANDELS; [Grogin et al.[[20Il) l, mid^nfrared 
{GOODS-Spitz er; PI: M. Dickinson; se e MagnelH et al. 2011), 
(sub)mm (e.g., [Oliver et al.|20I2|[BOTys et al.|2003nWeifi et alT] 
2009) and radio ( [Morrison et al.|20I0[|Miller et al.|2008| l obser- 
vations. 

We present here the combination of the PEP and GOODS- 
H observations of the GOODS fields. This combination pro- 
vides the deepest observations of the GOODS fields obtained by 
PACS. In particular, in the GOODS-South field the combination 
of these observations is not limited by the exposure time but by 
confusion. In this field, we thus obtain the deepest blank field 
observations achievable with PACS onboard the Herschel space 
observatory. In this paper we present in detail the data analysis 
method used to produce the publicly available PEP/GOODS-H 
maps and catalogue^ Then we use these deep observations to 
constrain the PACS-lOOyum and -160 ;um confusion noises and 
number counts, and to study the evolution of the infrared lu- 
minosity function and of the star-formation rate history of the 
Universe up to z ~ 2. 

The paper is structured as follows. Observations are de- 
scribed in Section |2] Section [3] presents the method used to 
produce the PACS maps. In Section]?] we present our source 
extraction methods and contents of the released package. In Sec- 
tion|5] we estimate the PACS-100 and -160 jum confusion noise. 
Number counts are presented in Section|6]while in Section|7]we 
derive the infrared luminosity functions of galaxies as well as the 
star-formation rate history of the Universe up to z~2. Finally, 
we summarise our results in Section [8] Throughout the paper 
we use a cosmology with Hq-II km s"' Mpc"', fl^ = 0.73 and 
Qm = 0.27. 



2. Observations 

The PACS maps and catalogues used and released in this paper 
are obtained from the combination of the PEP and GOODS-H 
observations of the GOODS-N and GOODS-S fields (see Table 
[TJ. Both programmes have observed the GOODS fields using 
the scan mode of the PACS photometer on board Herschel. This 
mode consists of slewing the spacecraft back and forth along 
parallel lines at a constant speed of 20" sec"'. Using this scan 
mode, astronomical observing requests (AORs) were designed 
to observe the GOODS fields in both nominal and orthogonal 
directions. To reach the desired depth many AORs per field were 
required. The central position of each AOR was dithered by ~ 8" 
in order to improve the spatial redundancy of the data. 

Using this strategy, the PEP and GOODS-H surveys cov- 
ered the entire 1 1' x 17' GOODS-N field with PACS at 100 and 
160 //nj^ The total observing time (i.e., including overheads) 
in GOODS-N was 25.8 and 124 hours for PEP and GOODS-H, 
respectively (Table[l]). 



^ http : //www . mpe .mpg . de/ir/Research/PEP/public_data_releases . p 

GOODS-H has also covered the entire GOODS-N field with 
SPIRE at 250, 3 50 and 500 yum SPIR E observations and catalogues 
are described in [Elbaz et al.[ ^2011^ and are publicly available at 
http://hedam.oamp. fr/GOODS-Herschel 
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Magnelli et al.: PEP/GOODS-H far infrared survey 
Table 1. Properties of the Herschel observations combined for tlie PEP/GOODS-H data release. 



Field 


Survey 


RA 


Dec 


Wavelengths 


Size 


Time 






[Degree 


, J2000] 




[arcmin] 


[h] 


GOODS-N 


PEP 


189.22862 


62.23867 


100, 160 


11' X 17' 


25.8 


GOODS-N 


GOODS-H 


189.22862 


62.23867 


100, 160 


ll'x 17' 


124.0 


GOODS-S 


PEP 


53.12654 


-27.80467 


70, 100, 160 


ll'x 17' 


226 


ECDFS" (GOODS-S) 


PEP 


53.10417 


-27.81389 


100, 160 


30' X 30' 


32.8 


GOODS-S 


GOODS-H 


53.12654 


-27.80467 


100, 160 


lO'x 10' 


206.3 


GOODS-S'' 


PV 


53.12654 


-27.80467 


100, 160 


11' X 17' 


7.9 



Notes. 

" Prior to be combined, observations of the ECDFS were trimmed to match the GOODS-S layout of the PEP observations. 

These observations were taken during the Herschel performance verification (PV) phase with preliminary instrument settings and thus were 
appropriately underweighted for the map creation (Section|3j. 



The entire 11' x 17' GOODS-S field was observed by PEP 
with PACS at 70, 100 and 160 yun]^ The total observing times 
were 113, 113 and 226 hours at 70, 100 and 160 //m, respec- 
tively, including 10 hours guaranteed time contributed by Her- 
schel mission scientist Martin Harwit. The Extended Chan- 
dra Deep Field South (ECDFS), containing in its centre the 
GOODS-S field, has also been observed by PEP with a total ob- 
serving time of 32.8 hours over a 30' x 30' region. Observations 
of ECDFS are included in our combined maps, trimmed to match 
the GOODS-S layout of the PEP observations. GOODS-H has 
observed a 10' x 10' region centred on the GOODS-S field with 
PACS at 100 and 160 ;um. The choice of covering only a fraction 
of the GOODS-S field was made in order to obtain, in a reason- 
able amount of time, a PACS- 100 fim map close to the confusion 
limit of Herschel. The GOODS-H observations of the GOODS- 
S field lasted for a total of 206.3 hours. Finally in our combined 
maps we also include 7.9 hours of observations of the GOODS-S 
field taken with preliminary instrument settings during the Her- 
schel performance verification phase (Table[T]i. 

Because PEP and GOODS-H observations were executed us- 
ing the same observing mode, a combination of these data sets 
can easily be performed. Of course, due to the different layout 
of PEP and GOODS-H observations in GOODS-S, this field is 
not homogeneously covered at 100 and 160jum. As illustrated 
by Fig.[T] a 10' x 10' region centred on the GOODS-S field has 
higher PACS -100 and 160yum coverage than the outskirts. In 
contrast, the GOODS-S 70 /im coverage is uniform across the 
field, with the exception of fall-off at the edges. In the rest of 
the paper, the centred region of the GOODS-S field with ultra- 
deep observations is referred as "GOODS-S-ultradeep", while 
the outskirts are referred to as "GOODS-S-deep". This com- 
bined data set provides the deepest Herschel far-infrared obser- 
vations of the GOODS-N and -S fields. In the rest of the paper 
this combined data set is referred to as the "PEP/GOODS-H" 
observations. 

Our combined PACS-lOO/im maps reach a total observing 
time per sky position of ~ 2.6 hours, ~ 2.6 hours and ~ 10.0 hours 
in GOODS-N, GOODS-S-deep and GOODS-S-ultradeep, re- 
spectively. The PACS- 160 //m maps reach a total observing time 
per sky position of ~ 2.6 hours, ~4.7 hours and ~ 12.1 hours in 
GOODS-N, GOODS-S-deep and GOODS-S-ultradeep, respec- 
tively. In the GOODS-S field, the PACS-70yum map reaches a 
total observing time per sky position of ~ 2.1 hours. 




s'as-oo' a'aa^ao' 3''32"oo' 

Right Ascension 
Center: R.A. 03 32 30.37 Deo -27 48 16-Bl 

Fig.l. PACS-lOOyumcoveragemapoftheGOODS-S field. Duetothe 
different layout of the PEP and GOODS-H observations, the 10' x 10' 
region centred on the GOODS-S field has much deeper PACS-lOO/jm 
observations than the outskirts. We observe the same pattern in the 
PACS-160;um coverage map of the GOODS-S field. 



We note that a detailed description of the observational 
strategies adopted by PEP and GOODS-H is provided in Lutz| 
et al. ( 201 1 1 and Elbaz et al. ( 201 1 1, respectively. 



' SPIRE-250/350/500yum observations o f the GOODS-S fi eld have 
been performed by the HerMES survey ( Oliver et al.||2012l l. These 
observations are publicly available through the Herschel Database in 
Marseille (HeDaM) at http : //hedam . oamp . f r/HerMES 



3. Map creation 

Observations were reduced using the standard PACS photometer 
pipeline (Wiep recht et al.|["2009| l and some custom procedures, 
all implemented within the HIPE environment in the Herschel 
common science system (HCSS). This data reduction process is 
described in more detail in |Lutz et al.| ( |20lT) and |Popesso et al.| 
( 2012) 1. Here we summarise the main steps. 

The data reduction process starts at the AOR level and is 
based on the scanmap script of the PACS photometer pipeline. 
First, the pipeline flags bad or saturated pixels, converts detec- 
tor signals from digital units to volts, finds pixels affected by 
short glitches and replaces their values using a standard interpo- 
lation method, and finally it applies a recentring correction to the 
pointing product of Herschel using reference positions of 24 yum 
sources with accurate astrometry (see |Lutz et al.|20lT] l from deep 
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Fig. 2. (Top left) PEP/GOODS-H map of the GOODS-north field (11' x 17') at 100 yUm. (Top right) PEP/GOODS-H map of the GOODS-south 
field (1 r X 17') at 100 yum. Contours correspond to exposure times greater than 0.5 (at the edge), 3, 6 and 9 (centre) hrs/pix. (Bottom left) Colour 
composite image of the GOODS-north field at 24pm (blue), 100 yum (green) and 160y[im (red). (Bottom right) Colour composite image of the 
GOODS-south field at 24 pm (blue), 100 pm (green) and 160 pm (red). The 24 pm images (PI: M. Dickinson) were obtained by the Spitzer Space 
Telescope, while the 100 and 160 /im images were obtained by the Herschel Space Obsen'atory. In the colour composite images, sources with 
contribution from the 24+lOOyum, 24-1-1 60 yum, 1004-1 60 yum and 24-1-100-1-1 60 yum bands would correspond to a cyan, magenta, yellow and white 
colours, respectively. The relatively noisy edges of the colour composite images have been trimmed. 



observations with the Multiband Imaging Photometer (MIPS; 
[Rieke et al. 20 04 ) onboard the Spitzer Space Telescope. Thanks 
to these recentring corrections, positions of PACS sources mea- 
sured in a "blind" catalogue (see Section 4.2 1 are offset by less 
than 0.2" and have a RMS difference of ~ 1" with respect to the 



position of their MIPS-24jum counterparts (see also Lutz et al. 
feoil) 



We note that the MIPS-24yum astro me try, and therefore 
our PACS maps, match the GOODS ACS version 2 coordinate 



systerrj^ Second, the pipeline removes from the timeline of each 
AOR the "1//" noise which is the main source of instrumental 
noise in PACS data. For deep cosmological surveys such as ours, 
the PACS pipeline does so by using a high-pass filtering method 
which subtracts, from each timeline, a version of the timeline 
filtered by a running box median of a given radius (expressed in 

« see the GOODS ACS data release at MAST: 
http: //archive . stsci . edu/prepds/goods/ 
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readouts, i.e., in numbers of points of the timeline). The presence 
of sources in the timeline affects the high-pass filtering method 
by artificially boosting this running box median, i.e., leading to 



the subtraction of part of the source flux ( Popesso et al.|2 012). 
Consequently, sources have to be "masked" from the timelines. 
Using results from Popesso et al. ( 2012| l, we choose a mask- 
ing strategy based on circular patches at prior positions. This 
method reduces the amount of flux loss due to the high-pass fil- 
ter and more importantly leads to flux losses which are indepen- 
dent of the PACS flux densities ( Popesso et al.|2012| l. We note 
that this improved masking strategy differs from that adopted in 
Lutz et al. (2011 1. Timelines are masked at the position of 24 



fim sources with 5 24 > 60 juJy using circular patches with radius 
of 4", 4" and 6" at 70, 100 and 160 fim, respectively. This 
strategy allows us to mask almost all PACS detections and cor- 
responds to the masking of about ~ 12% of the timeline. A small 
number of resolved sources (mainly in the GOODS-S field) are 
masked using larger patches which are adjusted visually. After 
being masked, timelines are high-pass filtered using a running 
box median with radius of 12 readouts (24"), 12 readouts (24") 
and 20 readouts (40") at 70, 100 and 160 fim, respectively. These 
radii are chosen based on results from Popesso et al. (2012) . Al- 
though we use an optimized high-pass filtering strategy, PACS 
flux densities still have to be corrected for flux losses. The cor- 
rections are provided in [Popesso et al.| ( [2012| l for our high-pass 
filter setting, i.e., 13%, 12% and 11% at 70, 100 and 160 //m, 
respectively (see Section|4|. 

After the removal of the "1//" noise, AORs are flat-fielded 
and flux calibrated. Then, the map of each AOR is created us- 
ing a "drizzle" method (Fruchter & Hook] '2002) implemented 
in the HCSS. Because all AOR maps are projected on the same 
world coordinate system, they are coadded into a final map us- 
ing, as weight, the eff'ective exposure time of each pixel, appro- 
priately underweighting the performance verification phase data 
which were taken with preliminary instrument settings. Taking 
advantage of the large number of AORs, we create uncertainty 
maps from the standard deviation of this weighted mean in each 
pixel. Although, a "drizzle" method is applied for the creation 
of each AOR map, the final map still contains some correlated 
noise. Thanks to the high redundancy of our data at each sky po- 
sition, we are able to calculate the mean correlated noise across 
the map and within our PSFs (see |Lutz et al.|20l"T| l. These mean 
correlation corrections were taken into account while measuring 
uncertainties from the uncertainty maps (i.e., uncertainties have 
to be corrected upward by a factor ~1.5; Section pT^) . 

Figure|2]shows the combined PEP/GOODS-H PACS-100 fim 
maps of the GOODS-N (top left panel) and GOODS-S (top right 
panel) fields. Three-colour composite images of the GOODS-N 
(bottom left panel) and GOODS-S (bottom right panel) fields at 
24- 100-1 60 jum are also shown in this figure. 



4. Source Extraction 

At the resolution of PACS most of the sources in our fields are 
point sources (i.e FWHM~4.7", 6.7" and 11" at 70, 100 and 
160 yum, respectively). Therefore, we use PSF fitting to derive 
their flux densities. Two catalogues are derived using com- 
plementary approaches. First, we construct a catalogue using 
as priors the source positions expected on the basis of a deep 
24 yum catalogue. This provides good deblending of neighbour- 
ing sources, but will miss a few Herschel sources that are not 
detected at 24 /im (i.e., < 1% and <4% of the PACS sources in 
the GOODS-N and GOODS-S fields, respectively). Second, we 



provide a "blind" catalogue using PSF fitting without positional 
priors. 

4. 1 . Prior source extraction 

Starting from the positions of IRAC-3.6;um sources (Infrared 
Array Camera; FWHM ~ 1 .6") from the GOODS Spitzer Legacy 
Program (PI; M. Dickinson), we ex tract sources in the MIPS- 
24//m maps (Magnelli et al. 201l| p and then use the 24yum- 
detected sources (i.e., with 5 24 > 3o" ~ 20//Jy) as priors for the 
source extraction in the PACS maps. The main advantage of 
this approach is that it deals with a large part of the blending is- 
sues encountered in dense fields and provides a straightforward 
association between IRAC, MIPS and PACS sources (Magnelli] 
|et al.lpOTT ). The disadvantage of this method is that we must 
assume that all sources present in our PACS images have al- 
ready been detected at 24 yum. Based on the relative depth of the 
mid- and far-infrared images of the GOODS fields, [Magdis et aL] 
( 20111 have investigated this assumption. They found that in the 
GOODS-H observations of the GOODS-S field, less than 2% 
of the PACS sources are missed in the MIPS-24jum catalogue. 
In our combined dataset, PACS observations of the GOODS-N 
field are shallower than those analysed in 'Magdis et al. ( 2011| l, 
while MIPS-24 yum observations are equivalently deep (i.e., 3cr ~ 
20 juJy in both GOODS-N and GOODS-S). Consequently, in the 
GOODS-N field, the fraction of PACS sources with no MIPS- 
24 jum counterparts should be significantly lower than 2%. This 
result is confirmed by cross-matching our GOODS-N "blind" 
PACS and MIPS-24 yum catalogues using a matching radius of 
4" and 6"at 100 and 160 /^m, respectively: less than 1% of the 
PACS "blind" sources are missed in the MIPS-24 jum catalogue. 
In GOODS-S, the situation is somehow different as in this field 
our PACS observations are deeper than those analysed in Magdis| 
eraLl([20ri). By cross-matching our GOODS-S "bhnd' PACS 



and MIPS-24 /im catalogues, we find that 4.1%, 4.2% and 3.5% 
of the PACS sources have no MIPS-24 jt/m counterparts within 
3", 4" and 6"at 70, 100 and 160 //m, respectively. These frac- 
tions should be treated as upper limits, given that most of these 
sources have faint PACS flux densities (i.e., 80% with S/N<5) 



and thus could be spurious (see Section 4.4 1. To understand the 



spectral properties of PACS sources with no MIPS-24 yum coun- 
terparts, we study the typical PACS-to-MIPS flux ratio observed 
in the GOODS-S field (Fig. |3]l. While the bulk of the PACS- 
100 (-160)//m population lies well within the parameter space 
reachable by our MIPS-24 /im catalogue (i.e., below the dashed 
lines of Fig. |3]l, there is, at faint flux densities, a slight trun- 
cation of the high-end of the dispersion of the PACS-to-MIPS 
flux ratio. This truncation will translate into faint PACS-100 
and -160yum sources with no MIPS -24 //m counterparts. In ad- 
dition, galaxies with similar spectral properties to Arp220 (i.e., 
with high PACS-to-MIPS flux ratio and strong silicate absorp- 
tions) will also be missed in MIPS-24 yum catalogues, especially 
at z~0.4 and z~ 1.3 where the ISyum and 9.4 jum silicate ab- 
sorption features are shifted into the MIPS-24yum passband (see 
also |Magdis et al. 201 1 1. The existence of PACS sources with no 
MIPS-24yL(m counterparts will naturally introduce slight incom- 
pleteness (i.e., < 1% and <4% in the GOODS-N and GOODS-S 
fields, respectively) in our prior catalogues. However, because 
this incompleteness mostly appears at faint flux densities (i.e., 
S/N < 5), it should not be much larger than that introduced by 



' The MIPS-24 /im catalogues used here are slightly different from 
those released by |Magnelli et al.|j2011^ . See Section [43] for more de- 
tails. 
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Fig. 3. PACS-to-24 yum flux density ratio as a function of the PACS flux density inferred from tlie "blind" source catalogue cross-matched with the 
MIPS-24yum catalogue. Shaded regions show the space density distribution of galaxies at 70yum {top left), lOOyum {top right) and 160 yum {bottom) 
in the GOODS-S field. In this field, the MIPS 24yum catalogue reaches a 2>cr limit of 20yuJy, while the PACS-70, 100 and 160 yum data reach 3cr 
limits of 0.9, 0.6 and 1.3 mjy (vertical dotted black lines). The parameter spaces reachable by our prior source extraction method using the MIPS 
catalogue are located below the dashed black lines: PACS sources located above these lines will not be detected at 24 ^m. Orange triangles, red 
circles, light blue squares and dark blue stars show the PACS-to-24 yum flux density ratios of galaxies with Ljr = 10"^ Lq as predicted at different 
redshifts usi ng the spectral ene rgy distributions (SEDs) of normal (MS), starburst (SB) galaxies (|Elbaz et al.|201I| l, IRAS22491 (Berta 2005) 
and Arp220 ( |Silva et al.|1998| l, respectively. The right point of each track corresponds to z = 0.5, while other points correspond to increasing 
redshifts, with intervals of Az = 0.5 (arrows indicate the path followed for increasing redshift). Predictions for galaxies with lower or higher 
infrared luminosities will correspond to a global vertical shift of these tracks towards lower or higher PACS flux densities, respectively. Most of 
the PACS population lies well within the limits of the parameter space reachable by our MIPS 24pm catalogues, i.e., below the dashed lines. 



prior-free source extraction methods at such low S/N. Indeed, at 
faint flux densities, there are only small discrepancies, in terms 
of number of sources, between our blind and prior catalogues 
(see Appendix [a). We note that we cannot use IRAC priors to re- 
duce this incompleteness because the high IRAC source density 
(i.e., 2-4 priors per PACS beam) would force each far-infrared 
source to be deblended into several unrealistic counterparts. 

From the expected positions of sources, we fit empirical 
PSFs (see Section |43| to our PACS maps. PACS flux densities 
are then defined as the intensity of the scaled PSFs. The simulta- 
neous fit of nearby sources optimizes the deblending of their flux 
densities. As for a standard aperture measurement, the photom- 
etry of each source is corrected to account for the finite size of 
our empirical PSFs, i.e., they do not include the power contained 
in the wings of the real PSFs (see Section [43] ). Finally, because 
of flux losses from the high-pass filtering, additional corrections 
are applied to our flux density measurements (see Section[3]l. 

Flux uncertainties are estimated on residual maps. These un- 
certainties are defined as the pixel dispersion, around a given 
source, of the residual map convolved with the appropriate PSF 
This method has the advantage of taking into account, at the 



same time, the rms of the map (including correlated noise and 
inhomogeneities in the exposure time) and the quality of our fit- 
ting procedure. Because our far-infrared observations have been 
designed to reach the confusion limit of Herschel, our flux un- 
certainties also contain a part of the confusion noise left in the 
residual maps (i.e., confusion noise, due to sources fainter than 
PACS detection limits). In this context of complex combina- 
tion of instrumental and confusion noise, we test the accuracy 
of our flux uncertainties using extensive Monte Carlo (MC) sim- 
ulations (see Section |4.4| i. We find that the flux uncertainties 
inferred from the residual maps and the photometric accuracies 
inferred from the MC simulations are in good agreement. 

4.2. Blind source extraction 

PACS flux densities are also extracted using a "blind" PSF-fitting 
analysis, i.e., without taking, as prior information, the expected 
positions of sources. This blind source extraction is performed 
with Star finder (D iolaiti et al. 2000 a b). This code is appro- 
priate for the extraction of the unrevolved PACS sources, be- 
cause it was especially designed to obtain high precision astrom- 
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etry and photometry of point sources in crowded fields. Ba- 
sically, Starfinder carries out source extraction using three 
steps: (i) detection of candidate point sources; (ii) verification of 
the likelihood of the candidate point sources; (iii) fits of empiri- 
cal PSFs to the map, at positions of the candidate point sources. 

Input parameters of Starfinder are fine-tuned using MC 
simulations (See Section [4!4l i. Starfinder performed its PSF- 
fitting analysis with the same empirical PSFs as those used for 
the prior source extraction (see Section 4.3 1. Flux densities in- 



ferred by Starfinder are corrected for the finite size of the em- 
pirical PSFs and for the effect of high-pass filtering (see Sections 
[3] and 4.3 1. From the MC simulations we observe that, despite 



the fine-tuning, Starfinder still tends to overestimate the flux 
density of sources at low S/N levels (about 10% at S/N ~ 3). This 
behaviour is corrected using factors derived from the MC simu- 
lations. 

Starfinder being exclusively designed for point source ex- 
traction, the presence of a few extended sources (~8 sources per 
field) in the PACS maps affects its detection process: extended 
sources are split into sub-components. To fix this problem, we 
run Sextractor on these extended sources using Kron ellipti- 
cal apertures. Extended sources are identified using an empirical 
method exploiting the isophotal area vs flux parameter space. 
In this parameter space, extended sources have large isophotal 
areas compared to their flux densities. PACS flux densities of 
the sub-components of the extended sources are erased from our 
blind source catalogue and replaced by flux densities inferred by 
SExtractor. We note that the accuracy of the Sextractor flux 
densities was verified using the Sextractor residual maps. 

Flux uncertainties are estimated with Starfinder using our 
PACS uncertainty maps (see Section [3]l and empirical PSFs. 
These flux uncertainties are corrected for the finite size of the 
PSFs and for the effect of high-pass filtering. In addition, be- 
cause our PACS uncertainty maps do not account for correlated 
noise, these flux uncertainties have to be corrected upwards us- 
ing mean correlation corrections (Section|3]l. We note that after 
correction, the flux uncertainties inferred by Starfinder are in 
good agreement with the photometric accuracies inferred from 
the MC simulations. 

4.3. PSF reconstruction 

Empirical PSFs used for the prior and blind source extraction 
are derived directly from the PACS maps using Starfinder. 
A number of isolated point-like sources present in the maps are 
stacked and then normalized to unit total flux. Because of the 
limited number of isolated point-like sources, the wings of these 
empirical PSFs have limited S/N ratio. Therefore, they have to 
be truncated to smaller radii (9.6", 7.2" and 12" at 70, 100 and 
160 /im, respectively) and then renormalized to a unit total flux. 

Because of their finite extent, the empirical PSFs do not in- 
clude all the power contained in the wings of the real PSFs. 
Consequently, flux densities inferred from these PSFs have to 
be corrected, as would be done for any standard aperture flux 
measurement. Aperture corrections are derived by comparing 
the empirical PSFs with the reference in-flight PSF for our ob- 
serving mode, obtained from observations of the asteroid Vesta. 
First, the Vesta observations are manipulated to account for vari- 
ations of the position angle (P.A.) between our AORs, i.e., the 
Vesta observations are rotated to the P.A. of each AOR and then 
stacked. Second, Vesta observations are smoothed to match the 
slightly broader FWHM of the empirical PSFs (broader by a fac- 
tor 1.04, 1.02 and 1.01 at 70, 100 and 160 yum, respectively). 
These broader FWHM can be explained by residual pointing un- 



certainties left by our recentring procedure (see Section |3]l. Fi- 
nally, aperture corrections are defined as the fraction of the to- 
tal energy contained in these manipulated Vesta PSFs within the 
radii of the empirical PSFs. At 70, 100 and 160 yum, the em- 
pirical PSFs contain 77, 68 and 69 % of the total energy of the 
manipulated Vesta PSFs, respectively. 

4.4. Monte Carlo simulations 

The PEP/GOODS-H observations have been designed to reach 
the confusion limit of the Herschel Space Observatory at 100 
and 160/im. Flux uncertainties are therefore a complex combi- 
nation of instrumental and confusion noise. In order to estimate 
these complex flux uncertainties and to characterize the quality 
of our catalogues we perform extensive MC simulations. 

We add simulated sources to our PACS 70, 100 and 
160yL(m images with a flux distribution, approximately match- 
ing the measured number counts (Berta et al.||2010| |2011| l. 
The flux densities of the faintest simulated sources are de- 
fined as the PACS flux densities (i.e., S a) expected for the 
faintest MIPS-24Aim (i.e., ^24 = 20AzJy) sources, i.e., 5™" = 
min(5'24) X mean(5',i/5'24), where mean(5'^/5' 24) is 10, 20 and 30 
ai A- 70, 100 and 160 //m, respectively (see Fig.[3]l. To preserve 
the original statistics of the image (especially its source num- 
ber density), the number of simulated objects added at one time 
is kept small (i.e., 20 sources). To recreate the clustering prop- 
erties of PACS sources, simulated objects are positioned, with 
respect to the MIPS-24yum sources, to reproduce the distance to 
the closest neighbour distribution observed in the MIPS -24 /im 
prior catalogues. Simulated sources are created using the manip- 
ulated Vesta PSFs which contain the wing of the real PSFs (see 
Section |43] l. On these simulated images, we then perform both 
our blind and prior source extraction using the empirical PSFs 
and compare the resulting flux densities to the input values. To 
improve the statistics, this process is repeated a large number of 
times using each time different positions and fluxes for the sim- 
ulated sources. For each field and at each wavelength, a total of 
20,000 artificial sources are used. Figure |4] shows, as an exam- 
ple, results from the MC simulations performed in the GOODS 
fields using the prior source extraction method. In both GOODS 
fields, these MC simulations correspond to the central deepest 
regions (i.e., avoiding the rather noisy edges of the GOODS- 
N field and concentrating on the GOODS-S-ultradeep part of 
the GOODS-S field). MC simulations on the outskirts of the 
GOODS-S field (i.e., GOODS-S-deep) have been run, but are 
not shown here. In any case, results from these MC simulations 
are very similar to those for the GOODS-N field, as expected 
from the similarity between the GOODS-S-deep and GOODS-N 
average exposure times. 

From these MC simulations we derive three important quan- 
tities: the completeness; the contamination; and the photometric 
accuracy of our catalogues as a function of flux density. Com- 
pleteness is defined as the fraction of simulated sources extracted 
with a flux accuracy better than 50%. The contamination is de- 
fined as the fraction of simulated sources introduced with S < la- 
hut extracted with S > 3cr. The photometric accuracy is defined 
as the standard deviation of the (5 out - 5in)/5out distribution as 
a function of 5 out (blue lines in Fig. [4]). These photometric accu- 
racies have the advantage of taking into account simultaneously 
nearly all sources of noise, i.e., including confusion. Table |2] 
summarises these quantities for both the blind and prior source 
extraction approaches. 

From these MC simulations we conclude that the prior and 
blind PSF-fitting methods perform accurate extraction of PACS 
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Fig. 4. Results of MC simulations in the GOODS-North and -South fields using the prior source extraction method. These MC simulations 
correspond to the central deepest regions of each field, i.e., avoiding the rather noisy edges of the GOODS-N field and concentrating on the 
GOODS-S-ultradeep part of the GOODS-S field. Blue lines represent the average photometric accuracy defined as the standard deviation of the 
(Soai - S;„)/S„^, distribution in each flux bin (after 3cr clipping). Red lines show the mean value of the (S^ui - Si„)/S„ut distribution in each flux 
bin. Inset plots show the fraction of artificial sources detected in the image (i.e., completeness) as a function of input flux (orange plain histogram) 
and the fraction of spurious sources (i.e., contamination) as a function of flux density (striped black histogram). 



sources. The prior and blind source catalogues are character- 
ized by high completeness and low contamination levels, as well 
as by good photometric accuracy (i.e., <33%). We observe in 
Fig. |4] that, using our flux uncertainties, a S/N > 3 cut almost 
always translates into a photometric accuracy better than 33% 
(i.e., as expected from sources with S/N > 3). We can thus con- 
clude that our flux uncertainties are fairly accurate. However, 
we also note that in the GOODS-S field the photometric accu- 



racy is worse than 33% for flux densities below 0.6 mJy and 
2.0 mJy at 100 and 160 jum, respectively. At such faint flux den- 
sities our PACS 100 and 160 yum maps are likely to be aff'ected by 
confusion noise (i.e., cr^(100//m) ~ 0.15 mJy and 0-^(160 fim) ~ 
0.68 mJy; see Sect. |5]l that is not fully accounted in flux uncer- 
tainty estimated by our source extraction methods. This conclu- 
sion is further confirmed by the analysis of the (5 in - S oud/Cs 
distribution (os being the flux uncertainty inferred by our source 
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Table 2. Statistics of released catalogues, for regions above the specified exposure time. 



Field 
& bands 


f > 

exposure 

hrs/pixel 


5(10-)" 
mjy 


Num.'' 
S/N>3 


Num.'' 
S/N>5 


Compl.'' 
% at 3cr 


Contam.'^ 
% at 3cr 


Compl.'' 
% at 5cr 


Contam.'' 
% at 5o- 


5(80%y 
mJy 


Prior source extraction 


GOODS-N 100 


1.3 


0.32 


914 


574 


36 


14 


72 


5 


1.97 


GOODS-N 160 


1.3 


0.70 


838 


546 


24 


32 


62 


5 


4.56 


GOODS-S 70 


1.0 


0.30 


476 


245 


21 


24 


77 


2 


1.39 


GOODS-S 100 ultradeep 


6.4 


0.18 


472 


346 


21 


27 


65 


5 


1.22 


GOODS-S 160 ultradeep 


8.3 


0.43 


445 


313 


18 


38 


45 


14 


3.63 


Blind source extraction 


GOODS-N 100 


1.3 


0.32 


865 


596 


13 


37 


56 


7 


2.04 


GOODS-N 160 


1.3 


0.70 


867 


521 


16 


48 


53 


17 


5.88 


GOODS-S 70 


1.0 


0.33 


396 


205 


19 


19 


72 


11 


1.51 


GOODS-S 100 ultradeep 


6.4 


0.18 


513 


377 


21 


42 


74 


9 


1.02 


GOODS-S 160 ultradeep 


8.3 


0.43 


485 


368 


15 


47 


45 


8 


4.89 



^ Numbers listed in this table are only suitable for sources situated in the regions of the field with exposure time higher than that reported in 
this column.To convert the PACS coverage maps provided in the released package into hrs/pixel, users should multiply them by 1.06 x 10"^. 
The \a flux density levels have been computed from 10 000 random extractions on residual maps. 
" Number of sources above a given S/N threshold. 

'' Completeness is defined as the fraction of simulated sources with 5 in ~3(5)(T and extracted with a flux accuracy better than 50%, i.e., 

-0.5<(5e„,-5i„)/5„„,<0.5. 
" Contamination is defined as the fraction of simulated sources introduced with 5in < 2cr but extracted with 5 out > Scr. 
^ Flux densities above which our catalogues are 80% complete. 



extraction methods). Indeed, in all but the GOODS-S 100 and 
160yum fields, the {S^ - Soudjcr^ distribution follows, as ex- 
pected, a Gaussian distribution with a dispersion almost equal to 
1. Instead, in the GOODS-S 100 and 160 yum fields, we find that 
the (5 in - Soudlcs distiibution follow a Gaussian distribution 
with a dispersion of ~ 1 .2, indicating that our flux uncertainties 
are underestimated and do not fully account for confusion noise. 
Therefore, we recommend caution in using flux densities lower 
than 0.6 mJy and 2.0 mJy at 100 and 160 ^m, respectively. 

From this analysis we conclude that in the GOODS-N field 
the combined PEP/GOODS-H data are 3.0 (1.16) and 2.8 (1.16) 
times deeper than the PEP (GOODS -H) data at 100 and 160 //m, 
respectively. Similarly, in the GOODS-S-ultradeep subarea cov- 
ered by both projects, the combined PEP/GOODS-H data are 2.3 
(1.5) and 2.0 (1.5) times deeper than the PEP (GOODS-H) data 
at 100 and 160jum, respectively. 

4.5. Content of the PEP/GOODS-H released package 

The PEP/GOODS-H released packag^ contains the scientific, 
uncertainty and coverage PACS maps. We reiterate that the 
GOODS-S 100/160 /im coverages and noise levels are inhomo- 
geneous due to the combination of observations with different 
layouts (see Section[2]i. In contrast, the coverages and noise lev- 
els for GOODS-N 100/160 //m and GOODS-S 70;um are fairly 
homogeneous, except for degradation near the edges. 

The released package contains PACS blind and prior source 
catalogues, down to 3cr significance. Because our source ex- 
traction methods might be inaccurate on the noisy edge of the 
PACS maps, we crop these regions from our catalogues. In or- 
der to track coverage variations across the fields, we provide 
for each source its exposure time (or equivalent) in each pass- 
band. The completeness and contamination levels inferred from 
the MC simulations are part of the released package. To use 
them, users should restrict their sample to sources with exposure 
time greater than that quoted in the second column of Table [2] 



Finally, we remind users that because our flux uncertainties do 
not fully account for confusion noise, sources with flux densi- 
ties below 0.6 m/y at 100 /um and 2.0 m/y at 160 jim have to be 
treated with caution. 

Calibration factors used to generate the final PACS maps are 
derived assuming an in-band SED of vSy = constant (see the 
Herschel data handbook for more details). Thus, some moder- 
ated colour corrections (i.e., < 7%) might need to be applied 
to our flux density measurements, as in-band SEDs of distant 
galaxies could be different from those assumed here. However, 
because these colour corrections depend on the redshift of the 
source, we decided not to apply any colour correction to the re- 
leased catalogues. 

PACS blind catalogues have been cross-matched to our 
MIPS-24//m catalogues using a maximum likelihood analysis 
( Ciliegi etal.|200 1 [[Sutherland & Saunders|I992! ). This method 
simultaneously accounts for fluxes and positions of MIPS -24 yum 
counterparts as well as positional errors in both the PACS and 
MIPS samples. The cross-identification of PACS and MIPS 
sources is included in the released package. In appendix |A] we 
compare the PACS blind and prior source catalogues using this 
cross-identification. 

Following the results of Elbaz et al. ( |201 1 1 and Hwang et al 
( [2010) 1, we provide for each source of the blind and prior cata- 
logues, its "clean index". Because this "clean index" is a mea- 
sure of the number of bright neighbours for a given source in all 
passbands, it supplies information on the potential contamina- 
tion of its flux densities. Here, "bright" neighbours are defined 
as sources brighter than half of the flux density of the source 
of interest, and closer than 20", 6.7" and 11" at 24, 100 and 
160, respectively. Having counted the number of neighbours of 
a given source (i.e., Neib24, Neibl88 and Neibl6Q), its "clean 
index" is given by 



clean_index = Neib24 + NeiblSO x 10 -i- Neibl60 x 100, 



http : //www.mpe .mpg . de/ir/Research/PEP/public_data_releases . php 



(1) 
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We note that the radius of 20" at 24 fim was defined in order to 
provide a "clean region" to the SPIRE-250;um galaxy popula- 
tion (FWHM- 18"; seeElbaz et al.||2011|). Therefore a "clean 
index" corresponding to Neib24 < 1 would provide a very con- 
servative selection for accurate PACS flux density estimates. 
Thus, we recommend the use of this conservative criterion (i.e., 
Neib24 < 1; 36% and 33% of the PACS sources in the GOODS- 
N and GOODS-S fields, respectively) not to select PACS sources 
with accurate flux densities but only to test whether any scien- 
tific results obtained with the full catalogue do not change when 
this criterion is applied. 

We stress that the MIPS-24yum catalogues released here 
are slightly different from those released by [Magnelli et al.| 
^MJ\ . These new MIPS-24//m catalogues are generated us- 
ing, as prior information, newer versions of the IRAC catalogues 
which are publicly available as part of the GOODS-H data re- 
lease. These MIPS-24jum catalog ues are identical to those used 
for the GOODS-H data releas^( Elbaz et al.|201 1] ). 

Finally, we note that the PEP/GOODS-H released package 
also contains ancillary products, i.e., the uncertainty and cov- 
erage maps, the empirical PSFs used by our source extraction 
methods, the Vesta PSFs used to measure aperture corrections, 
the results of the MC simulations and the PACS residuals maps. 



5. Confusion noise 



Following the formaUsm of| Dole et al] ( p003] l, deep PACS ob- 
servations might be affected by two different types of confusion 
due to extragalactic sources. First, the photometric confusion 
noise due to sources below the detection limit 5 Mm which pro- 
duce signal fluctuations (i.e., cTc) within the beam of the PACS 
sources, i.e., 5ii,n/cr(, =q with ^ = 3 or 5. Second, the density 
confusion due to the high density of sources above Sum which 
increases the probability to miss objects that are blended with 
bright neighbours. Because our GOODS-S 100 and 160 //m 
maps are the deepest blank field observations obtained by the 
Herschel Space Observatory, they are the most suitable datasets 
to estimate the PACS-100 and -160 jum photometric and density 
confusion noise. These estimates are presented in this section 
while the PACS-70jum confusion noise is derived and discussed 
in |BertaetaL| ( |20n] fliey find that the GOODS-S PACS-70//m 
observations are too shallow to be affected by either type of con- 
fusion). 

The photometric confusion noise is estimated empirically 
following the procedure described in [Frayer et al. ( 2006 ) and 
|BertaetaL] ( |2011[ ). First, we built several GOODS-S maps using 
only a fraction of the Herschel observations available. These 
partial-depth maps regularly probe the logarithmic exposure 
time parameter space from ~ 0. 1 hrs/pixel to ~ 10 hrs/pixel. Sec- 
ondly, we performed prior source extractions on these partial- 
depth maps and produced the corresponding residual maps re- 
moving sources with S >5']iin, with 5'iiin/o"c = ^ using ^ = 5. 
Thirdly, we measured the total noise ct-y in these residual maps 



using the procedure described in Sect. 4.1 i.e., ct-y is defined as 
the pixel dispersion of the residual map convolved with the Vesta 
PSF. Finally, we estimated cr,. by analyzing the variation of ct-y 
as a function of exposure time (see Fig. |5]l. This procedure was 
iterated until convergence at S'lim/cc = 5 was reached. 

Noise in partial-depth maps with short exposure time is dom- 
inated by the instrumental noise o"/ and therefore decreases as 
r^-^ . In contrast, noise in partial-depth maps with long exposure 
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Fig. 5. Noise in tlie PACS-100 /jm {top panel) and PACS-160ytim 
(bottom panel) maps as a function of exposure time. Empty circles 
represent the total noise ctj of the maps. The total noise is fitted with 
two noise components added in quadrature: an instrumental noise ctj 
component following a r"-^ trend (dashed line and empty squares) and 
an constant confusion noise component. The dotted lines present 
the two components fitted to ctt. The dot-dashed lines present the two 
components fitted to ctt made in [Berta et aT] (j20TTJ and illustrate the 
smaller exposure time range probed in their study. 



time departs from the f"**'^ trend and is a combination of instru- 
mental noise and confusion noise. Assuming that both compo- 
nents follow a Gaussian distribution, ct-y can then be approxi- 
mated by cTj - ^cr^ + cr^. Because cr^. does not depend on the 

exposure time, one can thus estimate cr^ by fitting ctj with a two 
components function, i.e., cr^ and crj oc f 

We find a photometric confusion noise cr^ of 0.15 and 
0.68 mJy in the PACS-100 and -160 //m passbands, respectively. 
The total noise of the PEP/GOODS-H observations is thus nearly 



http : //hedam. oamp . fr/GOODS-Herschel 



dominated by the photometric confusion noise at 100 yum (cTj""' = 
0.1 1 mJy) and it is fully dominated by the photometric confusion 
noise at 160 jum (cTj'^" = 0.21 mJy). In both bands, this signifi- 
cant contribution of the confusion noise to the total noise affects 
our flux uncertainty estimates: our flux uncertainties (Icr ~ 0. 18 
and 0.43 mJy at 100 and 160 //m, respectively) do not fully ac- 
count for confusion noise and sources with flux densities below 
0.6 mJy at 100 /im and 2.0 mJy at 160 fj.m have to be treated with 
caution (see Sect. |4.4[). 
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Fig. 6. PACS 100 and 160 yum differential number counts, normalized to the Euclidean slope {dN/dS oc 5 "^). Filled and open symbols 
show flux density bins above and below the 80% completeness limit, respectively. Grey shaded areas present estimates obtained using pre- 
Herschel observations. Blue shaded areas present estimates obtained using PACS observat ions jBerta et al.|2011^. Lin es represent predictions 



from baclcwards or forwards evolutionary model s (Lagache et al.||2004[ [Rowan-Robinson |2009| |Valiante et al.||2009[ |Le Borgne et al 
Franceschini et al.||2010[ [Gruppioni et al.||20T0l [Lacey et al.||2010[ [Marsden et al.||2011[ R^imati & van der Werf||2011[ [Niemi et"^ 



Bethermin etal.|2012> . 
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2012 



We note that our photometric confusion noise estimates 
are lower than those of [Berta et al.| ( |201 Ij ; cr^ ~ 0.27 and 
0.92 mJy at 100 and 160 jum, respectively). These discrepancies 
likely come from the fact that the observations used in |Berta| 



et al. (2011 1 were not deep enough to be dominated by the 
photometric confusion noise and therefore led to more uncertain 
estimates. This limitation is illustrated by the dot-dashed lines 
in Fig. [5] w hich represent the range of exposure time probed in 
[Berta et aLlpOTT) . 



catalogues, the density of sources is thus much higher than the 
density confusion limit defined by Dole et al. (2003 i.e., 16.7 



beams/source) but lower than those defined using more realistic 
assumptions (i.e., 3.5 beams/source). We stress that while an an- 
alytical estimate of the confusion density limit is useful, it should 
be used with caution as it does not account for the fact that our 
ability to separate pairs of sources depends on their S/N. In con- 
trast, empirical estimates through MC simulations take this effect 
into account. 



The confusion due to the high density of bright sources is 
usually defined as the flux limit 5 Mm at which the source density 
coiTesponds to 16.7 beams/source, i.e., the density at which 10% 
of the sour ces are separated by less than 0.8 x FWHM and thus 



blended" ( |Dole et aL|2603| ; where Q. = l.UxFW HM^ is the 



area of a beam). Using the PACS number counts, Berta et al. 
( |201 1| ) found that this density criterion corresponds to 5iim ~ 2.0 
and ~8mJy at 100 and 160 yum, respectively. These limits are 
thus greater than the detection limits derived from our MC sim- 
ulations. This discrepancy can be explained by the very conser- 
vative approach adopted by Dole et al. (2003 ). The density limit 
of 16.7 beams/source of , Dole et al., ( ,2003j corresponds to 10% 
of blended sources (i.e., separated by less than 0.8 x FWHM). 
This requirement translates into a blending-completeness of 90% 
while sources can be accurately extracted at lower completeness: 
at the 3cr level a catalogue has a typical completeness value of 
40-60% (i.e., 60-40% of blended sources). Using this more re- 
alistic assumption (i.e., measuring the density for which 40% of 
the sources are separated by less than 0.8 x FWHM), we infer a 
confusion density limit of ~ 3.5 beams/source. In our GOODS- 
S-ultradeep catalogues, the density of sources coiTesponds to 8 
and 4 beams/source at 100 and 160 /vm, respectively. In these 



6. Number counts 

Using the blind catalogues (i.e., without applying any "clean in- 
dex" selection) we inferred the PACS 100 and 160//m differ- 
ential number counts down to an unprecedented depths (PACS- 
70jum differential number counts are derived in [Berta et al.| 
201 ip. For that purpose, we applied the method described in 
Berta et al. (2010 201 l| l, which accounts for the incompleteness 
and contamination of our catalogues using results from the MC 
simulations (see also Chary et al. 2004; Small et al. |1995| l. In this 
method, observations with different depths cannot be treated si- 
multaneously. Therefore we divided our PACS observations into 
two sub-samples: (i) a ultradeep sub-sample containing sources 
in the GOODS-S-ultradeep field; and (ii) a deep sub-sample 
containing sources in the GOODS-S-deep and GOODS-N fields 
(these two fields having the same depth). The ultradeep sub- 
sample covers an effective area of 47 arcmin^, while the deep 
sub-sample covers an effective area of 327 arcmin^. 

The PACS 100 and 160 //m differential number counts in- 
feiTed from these two sub-samples are presented in Fig. [6] and 
provided in Table [A.l"| These differential number counts are nor- 
malized to the Euclidean slope (i.e., dN/dS ocS^^'^), expected 
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for a uniform distribution of galaxies in Euclidean space. Error 
bars include Poisson statistics, flux calibration uncertainties, and 



photometric uncertainties, as described in jBerta et al. (2011 



Filled and open symbols show flux density bins above and be- 
low the 80% completeness limit, respectively. Our differential 
number counts are compared with various estimates from the 
literature based on Spitzer or Herschel observations (for more 
details see Bert a et al. 2011) . There is good agreement between 
all these estimates over the flux density range in common. How- 
ever, thanks to the use of deeper PACS observations, our dif- 
ferential number counts extend to fainter flux densities than any 
previous estimates. This extension of ~0.5dex and ~0.2dex 
at 100 and 160/vm, respectively, allows us to resolve into in- 
dividual galaxies an even larger fraction of the cosmic infrared 



background (CIB). Using Figure 10 and CIB estimates of |Berta 
i.e 



et al. 
and 



(2011 



12.61+^^^ and 13.63+;^-^3 nWm-^sr-' at 100 
tJTfim, respectively, by power-law fitting of PEP data), 
we find that in the GOODS-S-ultradeep field our PACS obser- 



vations resolve ~15^_^f.^% and 



■75!{5% °f *e CIB at 100 and 



160yum, respectively. Furthermore, from the wealth of ancillary 
data available for the GOODS-S and -N fields (see Section |7]i, 
we can also study the redshift distribution of the PACS sources. 
Faint PACS 100 //m sources (i.e., S'loo [mJy] < 1.5) have a me- 
dian redshift of z = 1 -37^0 52 (errors give the interquartile range), 
while brighter sources (i.e., S 100 [mJy] > 1 .5) have a median red- 
shift of z = O.85+Q33. Similai-ly, faint PACS 160 jum sources (i.e., 
5 160 [mJy] <2.5) have a median redshift of z= 1.22''q4j while 
brighter sources (i.e., S leo [mJy] > 2.5) have a median redshift 
of z = 0.94+°-^2. lOOyum (160/im), sources at 0<z<0.5, 
0.5<z< 1.0,1.0<z<2.0 and z>2.0 contribute 24 (17)%, 36 
(33)%, 26 (28)% and 14 (22)% of the CIB resolved by our PACS 
observations. 

The PACS 100 and 160 jum counts are finally compared to 
predictions from backwards or forwards evolutionary models 
(|Lagache et al.||2004[ [R owan-Robin s6n|[2009} [Valiante et al. 



2009[|Le Borgne et al.||2009t |Franceschini et al.||2010[ |Grup 



pioni et al.||2010| |Lacey et al.||2010[ [Marsden et al.||201 1| |Rah 



mati & van der Werf||201 1[ |Niemi et al.||2012[ [Befliermin et al. 
2012[ l. Although most of the models reproduce the observed 
PACS 100 and 160 /im number counts fairly well, some of them 
can be ruled out because of significant discrepancies with our 
estimates. In particular, we note that the models of [Marsden | 
^taLl ( |20n] i, |Valiante et aT] ( [2009] l and |Niemi et"ari ( |2012! ) car> 
not reproduce the steep faint-end slope of the observed 100 and 
160 /im counts. From this cursory comparison, it is clear that 
ultradeep PACS number counts allow for a better refinement of 
the models and thus for better constraints on the evolution of 
star-forming galaxies. 



7. The infrared luminosity function 

Deep Herschel observations give us the opportunity to deter- 
mine the infrared luminosity function (LF) of galaxies with 
an unprecedented accuracy. Indeed, far-infrared observations 
provide more accurate infrared luminosity estimates than mid- 
infrared observations from Spitzer. for sources with multiple far- 
infrared detections (i.e., ~70% of sources in the PEP/GOODS- 
H fields; see Fig. |7]i, the SED-shape-LiR degeneracy is broken 
and infrared luminosity estimates are only limited by photomet- 
ric uncertainties; for sources with only one far-infrared detec- 
tion, uncertainties on the monochromatic-to-LiR conversion are 
significantly reduced compared to those provided by single mid- 
infrared detection (see Fig. |7]). The importance of far-infrared 



0.5 
" 0.4 

QJ 
t3 

^ 0.3 

X. 

a 

N 0.2 

b 0.1 

0.0 
0.4 

I 0.3 

o 0.2 

CO 

•fa 0.1 



0.0 



MIPS-24/im 
PACS-70/im 
PACS-lOO/xm 
PACS-160/im 




0.0 



Fig. 7. {Top panel) Uncertainties in determining Lir from monochro- 
matic observations (i.e., MIPS-24^m, PACS-70;um, PACS- 100 /zm or 
PACS-160;um bands) and the Dale & Helou (20021 SED library. These 
uncertainties are derived by taking the standard deviation of the log(LiR) 
distribution provided when normalizing all jPale & He lou (20021 SED 
templates to the same monochromatic flux density (±20% to account 
for typical photometric uncertainties, i.e., S/N~5) at a given observed 
wavelength (i.e., this observed wavelength depends on the band and 
redshift of interest). {Bottom panel) Fraction of PACS sources de- 
tected (i.e., S/N > 3) in only one of our PACS passbands (i.e., not with a 
70+100 //m, 70+160 yum, 100+160 yum or 70+ 100+ 160 ;/m detections) 
as a function of redshift. 



data increases further in cases where the mid-infrared may be 
contaminated by the emission from active galaxy nuclei (AGN). 

Taking advantage of PACS and SPIRE far-infrared observa- 
tions for several multi-wavelength fields, [Gruppioni et al.| ( |2013) 
derived the infrared LF of galaxies. Thanks to the large area 
covered by their observations (~ 2.5 deg^), they were able to ro- 
bustly constrain the intermediate and bright-end part of the in- 
frared LF up to z ~ 2 and the bright-end of the infrared LF up 
to z ~ 4. Here, we extend such study down to unprecedented 
depths using deeper PACS observations (~ 2 deeper in term of 
flux density than those of Gruppioni et al.). From these deeper 
observations, we are able to better constrain the faint-end and 
intermediate part of the infrared LF up to z ~ 2 and therefore to 
obtain better constraints on its redshift evolution. 

The GOODS-N and -S fields benefit from an extensive multi- 
wavelength coverage necessary to obtain redshift information 
for the PACS sources. In the GOODS-N field, we use a z + 
bands selected PSF-matched c atalogue created for the PEP sur- 



ve} 



10 



(Berta et al. 2010 2011), with photometry in 16 bands 
and a collection of spectroscopic redshifts (mainly from Barger] 
let al.|2008| ). In the GOODS-S field, we use t he GOODS-MUSIC 
z+K bands selected PSF-matched catalogue ( [Grazian et al. 2006 



|Santini et al.||2009| , with photometry in 15 bands and a coUec 
tion of spectroscopic redshifts. These multi-wavelength cata- 
logues also include photometric redshift estimates computed us- 
ing all optical and near-infrared data available (see |Berta et al 
lUTT] [Grazian et aL|[2006} [Santini et alT][2009l l. The quality of 



these photometric redshifts was tested by comparing them with 
the redshifts of spectroscopically confirmed galaxies. The high 
quality of these photometric redshifts is characterized by a rela- 
tively small scatter in Az/(1 +z) of 0.04 and 0.06 in the GOODS- 



publicly available at http : //www . mpe .mpg . de/ir/Research/PEP/ 
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N and -S fields, respectively ( |Berta et aT]|2011t [Santini et"aL] 
|2009] l. 

Multi-wavelength catalogues were cross-matched with our 
MIPS-PACS catalogues using their IRAC positions and a match- 
ing radius of 0.8" (i.e., approximately the HWHM of the IRAC- 
3.6 /im observations). In case of multiple associations (~ 10%), 
we select the closest optical counterparts. In the GOODS-N and 
GOODS-S fields, the common area covered by these catalogues 
is 164arcmin^ and 184arcmin^, respectively. In these regions, 
97% and 96% of the PACS sources have a multi-wavelength 
counterpart. Among those sources, 64% and 61% have a spec- 
troscopic redshift in the GOODS-N and -S fields, respectively. 
The rest of the sources has photometric redshift estimates. 

The total infrared luminosities (8-1000 fim) of PACS sources 
with redshift estimates were inferred by fitting their far-infrared 
flux densities (i.e., 70, 100 and 160 ;um) with the SED template 
library of Dale & Helou (2002), i.e., leaving the normalization 
of each SED template as a free parameter. For sources with only 
one far-infrared detection, infrared luminosities were defined as 
the geometric mean across the range of infrared luminosities 
given by all SED templates. As shown on Fig.jTj even in this case 
of only one far-infrared detection, the uncertainties in the in- 
ferred infrared luminosities are small, i.e., better than ~ 0.2 dex. 
To use the MIPS-24 yum flux density of the PACS sources during 
our fitting procedure does not change our results. Indeed, the 
LJ'i^'-'^^'^'^^/LJj^'^^ distribution has a mean value of 1 and a dis- 
persion of 3%. Naturally, for sources with only one far-infrared 
detection, the addition of the MIPS-24 yum flux densities allow 
us to break the SED-shape-LiR degeneracy and thus to reduce 
uncertainties on our infrared luminosity estimates. However, be- 
cause the fraction of PACS sources with single far-infrared de- 
tection is low (i.e., ~30%; see Fig. |7]) and because the MIPS- 
24 yum flux density might be affected by emission from an AGN, 
we decided not to use the MIPS -24 /im flux densities to derive 
LiR. We note that using the SED template library of [Chary &| 
|Elbaz| ( |200T] l instead of that of D ale & Helou] ( |2502l l to derive 
LiR (again leaving the normalisation as a free parameter), has no 
impact on our results. Indeed, the L™/L^°' distribution has a 
mean value of 1 and a dispersion of 13%. Figure[8]illustrates the 
detection limits of our PACS samples in term of Lir as a function 
of redshift. 

Uncertainties in determining Lir from monochromatic ob- 
servations (i.e.. Fig. |7]i are inferred assuming that the [Dale "&| 
[Helou ( 2002 ) SED library reproduces both the full range of mod- 
els appropriated for star-forming galaxies and the correct dis- 
tribution of SEDs within this population. Because neither of 
these assumptions are necessarily true, the absolute values of 
these uncertainties should be taken with caution. However, even 
with these limitations, uncertainties derived here are fully con- 
sistent with those inferred in Elbaz et al. (2011 1 by analyzing 



the mid-to-far-infrared SEDs (i.e., based on Spitzer, PACS and 
SPIRE observations) of a large sample of star-forming galaxies 
at0<z<2. 

The infrared (IR) LFs were derived using the standard 
1/^raax method ( [Schmidt[[I968[ ). The comoving volume of a 
given source is defined as y^ax = Vz„„„ - ^z,,,;,, where Zmia is the 
lower limit of the redshift bin being used, and Zm-dx is the max- 
imum redshift at which the source could be seen given the flux 
density limits of our observations, with a maximum value cor- 
responding to the upper Umit of the redshift bin. Here z^ax was 
defined by redshifting the [Dale & Helou[ ( [200 2 ) template fitted to 
the far-infrared flux densities of our sources until it fell below the 
detection limits of our PACS observations, or until Zms.x is greater 
than the upper limit of our redshift bin. For each luminosity bin. 
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Fig. 8. Infrared luminosities as a function of redshift for Herschel 
sources situated in the GOODS-S-uItradeep field (red squares) and sit- 
uated in the GOODS-S-deep and GOODS-N fields (open circles). The 
blue line on a white background shows the Lir limits above which the 
GOODS-S-ultradeep sample could be considered as a unbiased sample 
of the star-forming galaxy population at this redshift. These limits are 
inferred in Fig. [9] and "steps" correspond to the redshift bins used for 
the LP analysis. For the GOODS-S-deep and GOODS-N fields, these 
limits are shifted by ~ 0.2 dex towards higher Lir. 



the LF is then given by 

1 1 
= > , 

AL ^ y„,x,i X Wi 



(2) 



where Vmax,; is the comoving volume over which the /th galaxy 
of the luminosity bin could be observed, AL is the size of the 
luminosity bin, and w, is the completeness correction factor of 
the /th galaxy. The value of w, is given by the MC simulations 
(Sect. [4.4| i and depends on the flux densities of each source: w, 
equals 1 for bright PACS sources and decreases at faint flux den- 
sities due to the incompleteness of our PACS catalogues. For 
sources with multiple PACS detections it is defined as the maxi- 
mum of the three passbands, i.e., w/ = max(w™, w''^'^, w'^"). Be- 
cause we limit our LFs to infrared luminosities where w, > 0.5, 
none of our results strongly depend on these corrections. 

The minimum infrared luminosities that can be probed by 
the IR LFs in each of our redshift bins (i.e., 0. 1 < z < 0.4, 0.4 < 
z<0.7, 0.7<z<1.0, 1.0<z<1.3, 1.3<z<1.8 and 1.8<z<2.3) 
depend on the depth of our observations: at a given infrared 
luminosity, a large fraction of the galaxies has to be observ- 
able (i.e., small completeness correction; w, >0.5) over at least 
half of our redshift bin (i.e., small Vnjax correction; Zm-dx > 
^un' + (^3"' - ^un'^ll'^y this definition, it is clear- that 

observations with different depths cannot be treated simultane- 
ously. Therefore, as for the number counts, we first divided our 
PACS sources into two sub-samples: (i) a ultradeep sub-sample 
from the GOODS-S-ultradeep field; and (ii) a deep sub-sample 
from the GOODS-S-deep and GOODS-N fields. Then, for each 
of our redshift bins, we estimated the lowest infrared luminosi- 
ties probed by these sub-samples. 

At a given redshift, the minimum infrared luminosity ob- 
servable by PACS depends on the dust colour temperature of 
galaxies: at a given Lir, galaxies with warmer dust have brighter 
PACS flux densities. Therefore, we estimated, for each point of 
the Ldust-iiR parameter space its detectability by our PACS ob- 
servations using the SED templates of [Dale & Helou[ ( [2002| l, i.e. 
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Fig. 9. Selection limits introduced in the T^iist-Lm parameter space 
by our deepest PACS observations, i.e., in GOODS-S-ultradeep. Red 
continuous lines are these selection limits at different redshifts using 
our PACS 5(T detection limits. Each line corresponds to the central red- 
shift of our redshift bins, i.e., z = 0.25 (further left), z = 0.55, z = 0.85, 
z = 1.15, z = 1.55 and z = 2.05 (further right). The shaded area shows 
the local r[|„si-LiR relation found by Chapman et al. (2003 1, linearly ex- 
trapolated to 10'^ Lq. Dot-dashed lines show the T^^st-Lm relation in- 
ferred by Symeonidis et al.|^201 3 1 using a sample of high-redshift (i.e., 
0.2 < z < 1.2) //e«c/ie/-detected galaxies. Dashed black lines show, for 
each redshift, the lowest infrared luminosities probed by our ultradeep 
PACS observations without any dust temperature biases and yet popu- 
lated by star-forming galaxies. 



a dust temperature was assigned to each Dale & Helou template 
in a manner that is consistent with procedures used in Chapman] 
et al.| ([2003 ) to derive the Lir - Tdust relation. Then, assuming 



that the local Tdust-i'iR correlation of |Chapman et al.|(|2003|l re - 
mains the same at high-redshift (see also'Chapin et al. 2009 j ), 
we defined the minimum Lir of our IR LFs as the minimum Lir 
observable by PACS without any Ldust biases and populated by 
star-forming galaxies (i.e., within the T^^st-LiR correlation). In 
this analysis we used the 5cr limits of our PACS observations 
(i.e., where w, > 0.5) and the central redshift of our redshift bins 
(i.e,z„ax > 4ir+(<r -^lir)/2)- Hguregpresents the results 
of this analysis for our ultradeep sub-sampfe^^ The minimum Lir 
observable by PACS strongly increases with increasing redshift, 
and sources with hotter dust temperatures can be detected down 
to fainter infrared luminosities (red lines of Fig. [9]l. Combined 
with the expected positions of galaxies in the Ldust-i'iR parame- 
ter space, we can define the minimum Lir of our IR LFs (dashed 
lines of Fig. [9]l. The same analysis was performed for our deep 
sub-sample. Results are very similar but systematically shifted 
by ~ 0.2 dex towards higher Lir. We note that the modest evolu- 
tion of the Tdust-LiR relation with r edshift (see dot-da shed lines 
of Fig.|m Symeonidis et al. |2013| Hwang et al. |2010| Magnelh 
et al. in prep.) only has a minor effect on these limits, i.e., they 
are shifted by 0.05 -0.1 dex towards higher luminosities. We 
also note that using the dust temperature of galaxies situated on 
the main sequence (MS) of the SFR-M» plane, we would derive 
similar limits to those obtained using the Ldust-i^iR relation. In- 
deed, MS galaxies have Jdust of 27±3, 28+3, 29+3, 30+4, 32+4 



" To be consistent with Chapman et al. (2003 ), the T^^^t-^R relation of 
[Symeonidis et al.H2013) is obtained from their /60//100-Lir relation 
converted using the /60//100 and assigned dust temperature of each 
Dale & Helou template. 



and 34 ± 5 K at z = 0.25, z = 0.55, z = 0.8 5, 1.15, z= 1.55 and 
z - 2.05 (Magnelli et al. in prep; see also Magdis et al.|2012) . 

Uncertainties in the IR LFs depend on the number of sources 
per luminosity bin, on the photometric redshift errors and on 
the infrared luminosity errors. For the plot, they were defined 
as the quadratic sum of the Poissonian errors (ocl/ VA^) and 
errors computed from MC simulations which account for both 
photometric redshift and infrared luminosity uncertainties. The 



methodology of these MC simulations is described in Magnelli 
|etaL] ( |2009l|20ri] . 



Figure 10 



represents the IR LFs derived in six redshift bins 
using our ultradeep and deep PACS observations (Tables [A!2] and 
|A.3|l. We fit the IR LFs with a double power-law function sim- 



ilar to tha t used to fit the loca l IR LF which we also plot for 



reference ( [Sanders et al. 
and 



2003 



L for log(L/LQ) < Lk nee 

for log(L/Lo) >Lknee)- this fitting procedure, 



4> 



oc L 

the normalization (i.e., 0knee) and the transition luminosity (i.e., 
Lkn&e) of the double power-law function are left as free parame- 
ters. The shaded areas of Fig. [10| present the solutions compat- 
ible with the data within Icr. The evolution of 0knee and Lknee 
with redshift is presented in the upper left panel of Fig. [T2| and 
given in Table [A. 4| 

We compare our IR LFs with estimates made by [Magnelli 



[et al.j ( [2009| [201 1| ) using deep MIPS -24//m ob servations and 
//eric/ieZ-based estimates from Gmppi oiii et al.[ f2013) that use 
shallower PACS observations covering a wider effective area. 
At z< 1.8, we find good agreement with the Spitzer analysis 
of Magnelli et al. (2009 201 1 1. In particular, we observe good 
agreement, within the uncertainties, between our low luminos- 
ity extrapolations (based on 0knee and Lknee and a faint-end slope 
fixed at its z ~ value) and direct constraints obtained by Mag^ 



[nelli et aL ( 2009 2011 empty circles of Fig. 10 1. This agreement 
shows that even though our deepest PACS data do not allow us 
to probe luminosities far below the "knee" of the IR LFs, we 
obtain accurate low luminosity extrapolations, at least down to 
the infrared luminosities probed by Spitzer. The most signifi- 
cant difference between the present results and those of Mag^ 
'nelli et al.f (2009, '201 1) is observed at z ~ 2 and at the high- 



est infrared luminosities (Lir > 10'^ Lq). There, our new IR LF 
has a higher normalization. Because, the MIPS-24//m band can 
be affected by a significant contribution from an AGN, Mag-[ 
[nelli et al.[ ( |2009[ [201 l| l excluded from their sample all X-ray 
AGNs, i.e., sources detected in the GOODS-S/N Chandra cata- 
logues ( [Alexander et al.|r2003t [Lehmer et al.,,2005) with either 
Lx[0.5-8.0keV] > 3 x l(r^"erg s"' or a hardness ratio greater 
than 0.8 (I Bauer et al.[[2004| ). We excluded in the same manner 
all X-ray AGNs from our sample and found that these exclusions 
could not reconcile our IR LFs at z ~ 2. The observed difference 
is thus likely due to the fact that at z ~ 2, estimates from ,Mag^ 
[nelli et aT] ( [2009| 2011| l are affected by large pTe-Herschel un- 
certainties on the 24 yum-to-LiR conversion factors ( [Nordon et al.[ 
[2010, Elbaz et al. 2010 201 1 ; No rdon etal.|2012[ and reference 
therein). At this redshift, PACS estimates are more accurate. 

At all redshifts we observe very good agreement between 
our estimates and those of Gruppioni et al.l(2013). At high lu- 
minosities this agreement is very encouraging, because the data 
analysed by [Gruppioni et al.[ ( [2013| l sample a larger volume at 
shallower flux limits, and thus can more accurately measure the 
number densities for rare, bright sources. The only disagreement 
between our results and those of Gruppioni et al. (2013) appear 
at very low infrared luminosities, i.e., at luminosities not probed 
by both studies and which depend on the analytic function used 
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Fig. 10. Infrared luminosity functions estimated in six redsliift bins with tlie l/Vmax metliod. Red squares and black triangles show results 
from our ultradeep (i.e., GOODS-S-ultradeep) and deep (i.e., GOODS-S-deep and GOODS-N) samples, respectively. Dashed lines represent the 
best fits to these data points with a double power-law function with fixed slopes. The shaded areas span all the solutions of these fits which are 
compatible, within Icr, with our data points: the dark shaded parts of these areas highlight the luminosity ranges directly constrained by our PACS 
observations, while the light shaded parts highlight the luminosity ranges where our constraints rely upon extrapolations based on (^^nee and Li[„ee 
and a faint-end slope fixed at its z ~ value. Asterisks show the local reference, taken from |Sanders et aL| l2003 ), and the dotted line is the best 
fitted to these data points with our double power-law function with fixed slopes. Red dot-dashed lines are results from Magnelli e t al.|p(309| 
[201 1 ) using deep MIPS-24 yum observations. Magnelli et al. ( 2009 201 1 1 used the same double power-law function that is used here. To illustrate 
the infrared luminosity range constrained by Magnelli et al. (2009 2011 ), we show, as open circles, their lowest infrared luminosity bins. Blue 
triple-dot-dashed lines present results from |Gruppioni et al.,^20I3^ who used a difi'erent analytical function to fit their data points, in particular, a 
much shallower faint-end slope. 



to fit the data. |Gruppioni et al. ( 2013| l used a faint end slope 
shallower than that used here, i.e., -0.2 compared with -0.6. 
We find an evolution of (^knee (^knee = l0"2-5''*° '2 X (1 -H 

^)-I.5±0.7 for 2 < 1.0 and 0knee = 10-2"3±0.72 >^ + ^)-3.0±I.8 f^j. 
Z > 1.0) and Lknee (iknee = lO'^ ^^iO lO X (1 -H z)^ '^^" '' for Z < 1 .0 



and Lkne, 



10 



^knee " 
10.31±0.47 



X (1 -H zy 



4.2±1.2 



for z > 1 .0) in broad 



agreement with what was observed in [Magnelli et al.| ( [201 l[ l 
and Gruppioni et al. (2013 1. However, we note that in our study 
we find a stronger evolution of Lknee at z > 1 .0 than in [Magnelli] 
[etaL] ( |20n] i. 

With our ultradeep PACS data we can only constrain the IR 
LFs of galaxies down to Lir = 10" at z ~ 1 and Lir = 10'^ Lq 
at z~2. Consequently, at lower infrared luminosities, our 
constraints only rely on extrapolations based on (/>knee and Lknee 
and a faint-end slope fixed at its z ~ value. To test these extrap- 
olations, we take advantage of deep MIPS-24 yum observations 
using appropriate 24yum-to-LiR conversion factors. Thanks to 
recent Herschel studies, we know that the 24yL/m-to-LiR con- 



version factors depend not only on the MIPS -24 /im luminosity 
(as, e.g., in the Chary & Elbaz SED library) but also on the 
localization of galaxies with respect to the "main sequence" 
of star formation (MS; Iog(SFR)= a x Iog(M,) + C(z), where 
0.5 < Q- < 1.0; Elb az et aL|2011[|Nordon et al.|2012| l: galaxies 
situated on the MS have different 24yum-to-LiR conversion 
factors than galaxies situated above the MS ( [Elbaz e t al. 201 Ij 
jNordon et al.|[20l2) ). Consequently, we derive MlPS-based IR 
LFs using, on one hand, MS-based 24/im-to-LiR conversion 
factors and on the other hand, above-MS-based 24/vm-to-LiR 
conversion factors ( ]EIbaz et al.||201l"l Fig. [TT] i. The real IR 
LFs should be situated between these two estimates, but with a 
faint-end closer to the MS-based estimates (i.e., in a luminosity 
range dominated by MS galaxies) and a bright-end closer 
to the above-MS-based estimates (i.e., in a luminosity range 
dominated by above-MS galaxies). Up to z ~ 1.3, MS-based and 
above-MS-based estimates are in agreement: in this redshift 
range, MS-based and above-MS-based 24yum-to-LiR conversion 
factors are not significantly different. In contrast, at z > 1.3, 
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Fig. 11. Comparisons between the infrared luminosity functions derived using PACS observations and those derived using MIPS observations. 
PACS-based IR LFs are shown as dark/light shaded areas and dashed lines (see also Fig. \10\ . MlPS-based IR LPs are derived using MS-based 
24/jm-to-LiR conversion factors (black triangles) and above-MS-based 24yum-to-LiR conversion factors (red squares; see text for more details). 
The dotted lines show the local reference taken from [Sanders et ar| ( 2003l. Blue triple-dot-dashed lines present results from [Gruppioni et al.| 
( |2013) l. At low luminosities, the agreement between PACS-based and MlPS-based IR LFs confirms that PACS extrapolations towards lower 
infrared luminosities (i.e., light shaded areas) are reliable, at least down to the infrared luminosities probed by Spitzer. 



above-MS-based IR LFs have higher normalization than MS- 
based estimates: at z > 1.3, above-MS-based 24yL(m-to-LiR 
conversion factors are much larger than those for MS galaxies 
( |Elbaz et al.||20TT] [Nordon et aL][20T2 l). At z ~ 2 and at high 
infrared luminosities, the agreement found between PACS-based 
and above-MS-based estimates confirms the assumption that the 
bright-end of the IR LF is dominated by above-MS galaxies. 
At all redshifts and at low infrared luminosities (i.e., where 
MS galaxies should dominate the IR LFs), MS-based estimates 
are in perfect agreement with PACS extrapolations towards 
lower infrared luminosities (i.e., dashed lines and light shaded 
areas). We note that these PACS extrapolations also agree, 
within the error bars, with the above-MS-based estimates. All 
these agreements confirm that our PACS extrapolations are reli- 
able, at least down to the infrared luminosities probed by Spitzer. 

By integrating our PACS-based IR LFs we derive the evolu- 
tion of the comoving number density (top right panel of Fig. \T2\ 
and comoving infrar ed lii minosity density (IR LD, bottom panels 



of Fig 
iii 



12 



and Table 



A. 5 I of "faint" galaxies (i.e., 10^ L© < Lir < 
10"Lo)rLIRGs (i.e., 10" Lq < Ljr < 10'^ Lq) and ULIRGs 
(i.e., Lir > 10'^ Lq). Here, one has to keep in mind that even 
with our ultradeep PACS data we can only constrain the IR LFs 



of galaxies down to Lir = 10 Lq at z ~ 1 and Lir - 10 L© at 
z ~ 2, respectively. Therefore, number densities and IR LDs of 
galaxies below these limits rely upon extrapolations based on 
ipknee and L^nee and a faint-end slope fixed at its z ~ value. The 
redshift and luminosity ranges for which the inferred number 
densities and IR LDs rely upon these extrapolations are high- 
lighted by dashed lines and striped regions, respectively. Al- 
though these extrapolations seem to be corroborated by direct 
constraints from Spitzer (see Fig. Ill, we recommend caution 



when interpreting values not directly constrained by our PACS 
observations. We also emphasize that the LIRG and ULIRG des- 
ignations are used here strictly to segregate the luminosity bins, 
but not to imply physical properties. Indeed, Herschel studies 
have unambiguously revealed that high-redshift (U)LIRGs do 
not have the same properties as their local counterparts (e.g., EI 



baz et aI.|201ll|Wuyts et aI.l20TT| [Nordon et aI.|2012[rMagneUi 
et al.|2012| i. 



We find that the number density of ULIRGs strongly evolves 
from z = to z ~ 2, as it is multiplied by a factor ~ 800. Sim- 
ilarly, the number density of LIRGs significantly evolves with 
redshift: using direct constraints from PACS, we find that the 
LIRGs number density is multiplied by a factor ~ 25 between 
z = and z ~ 1, while relying on extrapolations we find that it is 
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Fig. 12. (Top left panel) Evolution of (pi^^^ and Lit„ee as a function of redshiift. {Top right panel) Evolution of the comoving number density 
of "faint" galaxies (i.e., 10^ Lq < Ljr < 10" U; black triangles), LIRGs (10" Lq < Lir < 10'^ Lq; orange diamonds) and ULIRGs (Lir > 10'^ Lq; 
red stars). Green circles show the number density of galaxies which are above the PACS detection limit of our ultradeep sub-sample (i.e., > L|r°, 
defined in Fig.|9j. The local reference is taken from Sanders et al. (2003). {Bottom left panel) Evolution of the total comoving IR energy density 
(black area) and the relative contribution of the "faint" galaxies (light blue area), LIRGs (orange area) and ULIRGs (red area). The black dot- 
dashed line shows the best fit of the total comoving IR energy density, i.e., defined by integrating the IR LPs best-fitting our data points (see 
dashed lines in Fig. |10| l. The green dashed line shows the best fit of the comoving IR energy density of galaxies which are above the PACS 
detection limit of our ultradeep sub-sample (i.e., > L"^'"). Black circles represent the total comoving IR energy density inferred by Gruppioni et al. 
(2013). The right axis of the panel shows the evolution of the obscured SFR density assuming that SFR [M© yr"']= 10"'" x Lir [Lq] for ajChabrier 
( 2003 1 IMF. {Bottom right panel) Evolution of the comoving IR energy density. Areas are the same as in the bottom left panel. The dotted lines 
represents the unobscured SFR density of the Universe (i.e., not coiTected for extinction; Cucciati et al.|2012| . The red dashed line on a white 
background shows the total SFR density of the Universe, defined as the sum of the obscured and unobscured SFR densities. Yellow diamonds are a 
compilation of SFR density estimates from Hopkins & Beacom ( 2006 ). Green circles and squares are the total SFR density of the Universe and the 
contribution of ULIRGs estima ted by|Murphy et al. (201 1) using deep MIPS-24 and -70 pm observations. Dark blue diamonds, stars and triangles 
are estimates from |Casey et al.| ( |2012^ for all galaxies, LIRGs and ULIRGs, respectively, using SPIRE (i.e., submillimetre) observations. One has 
to keep in mind that even with our ultradeep PACS data we only constrain the IR LFs of galaxies down to Lir = 10" Lq at z ~ I and Lir = 10'^ Lq 
at z ~ 2, respectively. Therefore, number densities and IR LDs of galaxies below these limits rely upon extrapolations based on 0ii„ee and Lk„ee and 
a faint-end slope fixed at its z ~ value. The redshift and luminosity ranges for which the inferred number densities and IR LDs rely upon these 
extrapolations are highlighted by dashed lines and striped regions, respectively. We note that these extrapolations are nevertheless corroborated by 
direct constraints based on deep MIPS-24 ;um observations (see Fig.[TTJ. 



multiplied by a factor ~ 30 between z = and z~2. Naturally, 
the redshift evolution of the LIRGs and ULIRGs number den- 
sities translate into the redshift evolution of their IR LDs. The 
IR LDs of LIRGs and ULIRGs are multiplied by a factor ~ 130 
(~ 1000) and ~ 40 (~ 45; based on extrapolations) between z = 
and z ~ 1 (2), respectively. The redshift evolution of the IR LD 
of ULIRGs is consistent with that found by |]VIurphy et al.| ( |20TT] l 
using deep MIPS-24/70;um observations (dark blue squares in 
the bottom right panel of Fig. 12 1. In contrast, at z ~ 2, .Magnelli 



[etaL] ( [20n] ) found an IR LD of ULIRGs lower by a factor ~ 3. 
This inconsistency reflects the discrepancies observed at z~2 
between their and our IR LF (see Fig. 



10 1 



Relying upon extrapolations to low infrared luminosities, we 
find that the total IR LD strongly increases from z = and z ~ 2, 
with most of the evolution happening before z ~ 1 .2 driven by the 
strong increase of the IR LD of LIRGs. At z ~ 1, LIRGs account 
for 50 + 26% of the total IR LD while ULIRGs contribute only 
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10 + 6% of it. In contrast, at z ~ 2 ULIRGs contribute 50 ± 24% 
while LlRGs contribute only 30 ± 20% of the total IR LD. 

Although the total IR LDs derived here are consistent, within 
the uncertainties, with those from Magnelli et al. (2011), there 
are some differences. At z~ 1, the total IR LD derived here is 
lower by a factor ~ 1.3 than that of [Magnelli et al. (201 1 1 while 
at z~2 it is higher by a factor ~ 1.3. More importantly, at z~ 
2, the ULIRGs contribution to the total IR LD found here (i.e.. 



50% ±24) is much higher than that found in Magnelli etal. (2011 



~ 17%). Naturally, this finding reflects the disagreement at z ~ 2 
between our two IR LFs. 

The total IR LDs inferred in'Gruppioni et al. ( 2013 ! black cir- 
cles in the bottom left panel of Fig. [T2| agree with our findings. 
Noticeably, their and our estimates are obtained by integrating 
IR LFs which significantly differ at infrared luminosities lower 
than those probed by Spitzer, i.e., [Gruppioni et al. (2013) used 
a much shallower faint-end slope (see Fig. [TTJ. This demon- 
strates that extrapolations towards very low infrared luminosi- 
ties do not aff'ect much the estimates of the total IR LDs. Indeed, 
unless the faint-end slope of the IR LFs is significantly steeper 
than that used here, most of the total comoving IR LD is emitted 
by galaxies with luminosities reachable with PACS (i.e., ~ 50%; 
see black dot-dashed and green dashed lines in the bottom left 
panel of Fig. 12 1 and Spitzer (i.e., ~ 75%). 

While the IR LDs derived here for all galaxies and for 
LlRGs are consistent with SPlRE-based results from |Casey| 
leTaLl (|20T2ji, the IR LDs of ULI RGs at z > 0. 8 deriv ed here 
are lower than those inferred by |Casey et al.| ( |2012[ ). This 
disagreement might be due to the large spectroscopic redshift 
incompleteness corrections applied to the IR LFs of 'Casey et al. 
(2012 ) at z>0.8 (i.e., corrections >50%). This seems to be 
confirmed by the good agreement, at z > 0.8, between the IR 
LFs of Gruppioni et al. ( 2013| i.e., also in agreement with our 



IR LFs) and the SPIRE-based IR LFs of Vaccari et al. (in prep). 
Finally we note that the total IR LD derived here, as well as the 
significant contribution of ULIRGs at z ~ 2, are consistent with 
results from [Murphy et al.| ( |201 \) . 

Assuming that the IR LD is totally produced by star forma- 
tion (i.e., without any contribution from AGNs), it can be con- 
verted into the obscured SFR density of the Universe using the 
relation of |K ennicuttl ( |1998[ ), scaled to a |Chabrier|p003| l initial 
mass functiori'*^ 



SFR[Moyr~i] 



10-'°xLir[Lo], 



(3) 



About 10-15% of our PACS sources are associated with X-ray 
AGNs that might significantly contribute to the IR LD derived 
here. However, Herschel studies of X-ray AGNs ([Rosario et"ar| 



2013[|Shao et al.|2010[pullaney et al.|2012[|Rosario et al.|2012 1 
have demonstrated that in the vast majority of cases (i.e., > 94%) 
the PACS flux densities are dominated by emission from the 
host galaxy and thus provide an uncontaminated view of their 
star-formation activities. Consequently, we assume that the IR 
LD derived here has no significant contribution from AGNs and 
can be converted into the obscured SFR density of the Universe 
(right axis of the bottom panels of Fig. [T2| ). Combined with 
the unobscured SFR density of the Universe derived by Cucciati 
et al. ( 2012| l using rest-frame UV observations, we can then infer 
the total SFR density and its evolution up to z ~ 2 (dashed red line 



on a white background in the bottom right panel of Fig. 12 1. The 
cosmic star-formation history strongly evolves from z = to z ~ 1 



and flattens at z > 1 . The unobscured SFR density accounts for 
about ~ 25%, ~ 12% and ~ 17% of the total SFR density of the 
Universe at z ~ 0, z ~ 1 and z ~ 2, respectively. The contribution 
of the unobscured SFR density and its evolution with redshift is 
consistent with the redshift evolution of the mean rest-frame UV 
dust attenuation (e.g., Cucciati et al.||2012 Tresse et al. |2007 



Finally, we note that the cosmic star-formation history derived 
here is fully consistent with the combination of indicators, either 



obscured or corrected for dust extinction, as compiled by Hop- 
|kins & Beacom|p006| l, [Seymour et al.| ( |2008| l and [Karim eTaT 
( 20 1 1 not shown in Fig.[12|l. 



Kennicutt ( 1998 i adopts a Salpeter initial mass function so that we 



divide his normalization by 1.72. 



8. Summary 

By combining observations of the GOODS fields from the PEP 
and GOOT)S>-Herschel key programmes, we obtain the deepest 
PACS far-infrared blank field extragalactic survey carried out by 
the Herschel Space Observatory. In particular, in the GOODS- 
S field the combination of these observations is not limited by 
the exposure time but by confusion. PACS flux densities are 
extracted from the maps using two complementary PSF-fitting 
approaches. Firstly, we extract PACS flux densities using, as 
prior information, the expected positions of the sources on the 
basis of deep 24-fim catalogues. Secondly, PACS flux densities 
are extracted "blindly", i.e., without positional priors. The ac- 
curacy of both approaches is tested through MC simulations. In 
the deepest parts of the GOODS-S field, these catalogues reach 
3cr depths of 0.9, 0.6, 1.3 mJy at 70, 100, 160 yum, respectively. 
From these catalogues we derive number counts down to these 
unprecedented depths, and determine the infrared luminosity 
function down to Ljr = 10" Lq at z ~ 1 and Lir - 10'^ Lq at z ~ 2, 
respectively. By integrating these infrared luminosity functions, 
we estimate the evolution of the SFR density of the Universe up 
to z ~ 2. We find that the obscured SFR density of the Universe 
strongly increases from z = to z ~ 1 and then increases more 
moderately up to z ~ 2.3. The obscured SFR density of the Uni- 
verse is dominated by the contribution of the LlRGs at z ~ 1 (i.e., 
50 ± 26%) while it is dominated by the contribution of ULIRGs 
(50 ± 24%) and LlRGs (30 ± 20%) at z ~ 2. 

Maps and source catalogues (> 3cr) are now publicly re- 
leased. Combined with the large wealth of multi-wavelength 
data available for the GOODS fields, they provide a powerful 
new tool for studying galaxy evolution over a broad range of 
redshifts. 
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Appendix A: Comparison between blind and prior 
catalogues 

Figure [AT| presents the comparison between the blind and prior 
catalogues in the GOODS-N and GOODS-S fields. These com- 
parisons are restricted to the region of the field with exposure 
time higher than those mentioned in the second column of Ta- 
ble |2] For the GOODS-S 100 and 160 //m maps, this restric- 
tion corresponds to the central deepest region of the field, i.e., 
GOODS-S-ultradeep. The cross-identification of PACS blind 
and prior sources is based on their MIPS-24/im positions, i.e., 
for the blind catalogues we use the PACS/MIPS-24/im cross- 
identification performed using a maximum likelihood analysis 
(see Section [43] l. For each field and passband we present the 
direct comparison of flux densities for sources in common to 
both catalogues; the distribution of unmatched sources in abso- 
lute number; and the distribution of unmatched sources in frac- 
tion relative to total in the given flux density bin. 

Flux densities extracted with these two independent source 
extraction methods are consistent with each other. At faint 
flux densities, there are however few outliers, with higher flux 
densities in the blind catalogues than in the prior catalogues. 
Such outliers are expected since blind source extraction is more 
severely affected by blending issues leading to the overestima- 
tion of PACS flux densities. 

In the GOODS-N field the absolute number of unmatched 
sources in the blind and prior source catalogues are very similar. 
At faint flux densities, the increase of the fraction of unmatched 
sources is likely explained by the increase of the incompleteness 
of our catalogues. Indeed, because our two source extraction 
methods are independent, at faint flux densities their incomplete 
samples might not fully overlap. 

The blind GOODS-S 100 and 160 //m source catalogues 
contain a larger number of unmatched sources than the prior 
source catalogues. Some sources in "excess" in the blind cat- 
alogues should correspond to sources effectively missed by our 
prior extraction due to the lack of MIPS-24 /im counterparts (see 
Section |4.1| i. However, a significant fraction of the sources in 
"excess" in the blind catalogue should correspond to spurious 
detections, as in this field the contamination of the blind cata- 
logues is supposed to be higher than that of the prior catalogues 
(see Table The GOODS-N 100 yum, GOODS-S 100 A^m and 
GOODS-S 160 jum blind catalogues also contain each one a 
bright source (i.e., > 8 mJy) missed by the prior extraction. Ex- 
amining these sources, it turns out that they are not missed be- 
cause of a lack of MIPS-24 /im priors, but likely correspond to 
"spurious" bright sources created by the blind extraction by over- 
deblending very bright and crowded regions of the field. Finally, 
we observe that the prior GOODS-S 70 /im source catalogue 
contains a larger number of unmatched sources as it reaches a 
slightly lower 3cr detection limit than the blind catalogue. 



Article number, page 19 of 



22 




m 

4 


GOODS-S 

U" ■ 

■ 




R_ ■ 


B ■ ■ 




Nb sources : 369 
1 1 



10 



e 



Prior (75) 
Blind (117) 



10 

mJy] 



S 10 

T3 

OJ 

■g 0.8 

c 0.6 

=) 

° 0,4 

c 

o 

0.2 
0.0 



Prior (75) 
Blind (117) 




10 
mJy] 



Li I 




Nb sources : 423 



nJy] 



s 



Prior (48) 
Blind (89) 




) [mJy] 



S 1.3 

u 

I 1.0 
-a 

0} 

c 0.6 

= 0.4 

c 

o 

u 0.2 
i- 

0.0 



J1 



1 



Prior (46) 
Blind (89) 



nJy] 




Nb sources : 337 



10 
[mJy] 



Prior (139) 
Blind (59) 



10 
[mJy] 



S 1.0 

-o 

OJ 

•g o.e 

c 0.6 
□ 

° 0.4 

c 

o 

0.0 



Prior (139) 
Blind (59) 




10 
[mJy] 



Fig. A.l. Comparison between blind and prior source catalogues. (Left panel) Direct comparison of flux densities for sources in common to both 
catalogues. (Centre panel) Distribution of unmatched sources in absolute number. (Right panel) Distribution of unmatched sources in fraction 
relative to total in the given flux density bin. 
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Magnelli et al.: PEP/GOODS-H far infrared survey 
Table A.l. PACS 100 and 160 yum number counts, normalized to the Euclidean slope. 



PACS lOOjum PACS 160jum 



GOODS-S ultradeep GOODS-N/S deep GOODS-S ultradeep GOODS-N/S deep 



^ centre 


Counts 


Err. 


Counts 


Err. 


c 

centre 


Counts 


Err. 


Counts 


Err. 


[mJy] 


[lOMeg 


-2mJy'-5] 


[lO^deg-'mJy'^] 


[mJy] 


[10" deg 


mJy'"*] 


[lOMeg 


-2 mJy'-^] 


0.56 


0.59 


0.10 


- 


- 


1.42 


3.59 


0.49 


- 


- 


0.71 


0.91 


0.14 


- 


- 


1.79 


4.91 


0.72 


- 


- 


0.89 


1.34 


0.20 


- 


- 


2.25 


6.49 


0.91 


4.92 


0.31 


1.13 


1.78 


0.28 


1.29 


0.10 


2.84 


7.56 


1.15 


6.48 


0.41 


1.42 


1.95 


0.36 


1.90 


0.13 


3.57 


8.62 


1.54 


8.80 


0.56 


1.79 


2.80 


0.51 


2.76 


0.19 


4.49 


9.89 


1.93 


11.85 


0.78 


2.25 


3.59 


0.70 


3.63 


0.26 


5.66 


15.08 


2.80 


14.36 


1.03 


2.83 


3.70 


0.85 


4.27 


0.33 


7.13 


15.37 


3.32 


15.55 


1.28 


3.57 


3.73 


1.03 


4.79 


0.42 


8.97 


13.30 


3.86 


16.68 


1.59 


4.49 


3.51 


1.23 


5.44 


0.55 


11.29 


18.79 


5.41 


18.72 


2.02 


5.66 


6.54 


1.93 


6.25 


0.72 


14.22 


14.87 


5.93 


17.15 


2.32 


7.13 


7.53 


2.48 


5.72 


0.84 


17.91 


10.63 


6.23 


16.86 


2.79 


8.97 


4.71 


2.54 


5.87 


1.02 


22.54 


24.76 


10.55 


19.48 


3.63 


11.29 


6.19 


3.18 


4.21 


1.04 


28.38 


21.99 


12.47 


17.99 


4.21 


14.22 


5.37 


3.65 


5.85 


1.42 


35.72 






15.81 


4.69 


17.90 






8.36 


1.99 


44.98 






11.07 


4.68 


22.54 


9.46 


6.44 


5.66 


2.05 


56.62 






11.82 


5.63 


28.38 






7.97 


2.65 













Table A.2. The infrared LF derived from the 1/ V,iiax analysis in the GOODS-S ultradeep field 



log(L;°/) - Xogd'X) log(0) 

[log(Lo)] [log(Mpc-3dex-l)] 


log(Li™) 

[log 


-log(Lr) 

(Lo)] [lo 


log(0) 

g(Mpc"^dex"')] 




0.1 <z<0.4 






0.4<z<0.7 




9.7- 
10.1 - 
10.5 - 
10.9- 


10.1 
10.5 
10.9 
11.3 


9 13 
^•■'■'-0.20 

_9 4^+0.13 

9 CQ+0.15 
^■-'^-0.25 

T 7Q+0.21 
■^•^^-3.29 


10.4 
10.8 
11.2 


- 10.8 

- 11.2 

- 11.6 


9 ac!+0.I3 
^"'-'-0.14 

-2 56+0-14 

^•-'"-0.15 
9 Q/c+0.16 




0.7<z<1.0 






1.0<z<1.3 




10.8 - 
11.2- 
11.6- 


11.2 
11.6 
12.0 


-2 55+013 
-3 05+°"' 

■^•^■^-0.32 


11.2 
11.6 
12.0 


- 11.6 

- 12.0 

- 12.4 


_9 74+0.17 
-0.18 

-3 09+° ''* 
-0.22 

_4 1 c+0.34 
^■^•'-4.15 




1.3<z<1.8 






1.8<z<2.3 




11.4- 
11.8 - 
12.2- 


11.8 
12.2 
12.6 


9 Ql+0.17 
^■^-"-O-IS 

-3 20+° '^ 
-3 96+0-25 

■^•^"-0.39 


11.9 
12.3 


- 12.3 

- 12.7 


T 1 c + 0.14 
•'■^-'-0.16 

-3 67+° '^ 

■'■"'-0.25 
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Table A.3. The infrared LF derived from the l/V,,,;,^ analysis in the GOODS-N/S deep fields 



log(Zjp^log(Z;f5 log(0) log(L;™)-log(Lff) log(0) 

[log(Lo)] [log(Mpc-3dex-l)] [log(Lo)] [logCMpc-^dex-')] 



0.1<z<0.4 0.4<z<0.7 



10.0- 
10.4- 
10.8 - 
11.2- 


10.4 
10.8 
11.2 
11.6 


_o 47+0.05 
^•^ ' -0.07 

-2 69+"°^ 

o 1Q+0.12 
-'■^^-0.18 

-4 09+0-29 


10.6- 
11.0- 
11.4- 
11.8 - 


11.0 
11.4 
11.8 
12.2 


-2 46+°-" 

9 0^+0.12 
^■^-'-0.12 

-'■^^-0.17 

-4 00+°-^° 

^■^ -0.31 




0.7<z<1.0 






1.0<z<1.3 




11.0- 
11.4- 
11.8 - 
12.2- 


11.4 
11.8 
12.2 
12.6 


7 CQ+0.10 

^•-'^-o.io 

_C! 14+0.11 

-2 OO+0.16 
-3.»»_Q20 

-4 83+''-3' 


11.4- 
11.8 - 
12.2- 
12.6- 


11.8 
12.2 
12.6 
13.0 


-2 98+"-'5 

^■^°~0.16 

T 'iQ+o.n 

-'-^°~0.18 

-4 47+0.25 
^■^'-0.40 

_4 QC+0.34 
^■^-'-4.95 




1.3<z<1.8 






1.8<z<2.3 




11.6- 
12.0- 
12.4- 


12.0 
12.4 
12.8 


o 1 Q+0.16 
■'•^^-0.15 

o 7Q+0.I6 
_4 ^4+0.20 


12.1 - 
12.5 - 
12.9- 


12.5 
12.9 
13.3 


-3 50+"-" 

-0.12 

-4 06+°-''' 

_4 QQ+0.25 
^■^^-0..54 



Table A.4. Parameter values of the infrared LF 



Redshift 


«i^ 




log(Lknee) 


log(0knee) 








[log(Lo)] 


[log(Mpc"3(jgx-l)] 


z~0 


-0.60 


-2.20 


10.48 + 0.02 


-2.52 + 0.03 


0.1 < z < 0.4 


-0.60 


-2.20 


10.84 + 0.06 


-2.85 + 0.04 


0.4 < z < 0.7 


-0.60 


-2.20 


11.28 + 0.12 


-2.82 + 0.14 


0.7 < z < 1.0 


-0.60 


-2.20 


11.53 ±0.15 


-2.96 + 0.18 


1.0 < z < 1.3 


-0.60 


-2.20 


11.71+0.14 


-3.01 ± 0.20 


1.3 <z< 1.8 


-0.60 


-2.20 


12.00 + 0.15 


-3.29 + 0.19 


1.8 <z<2.3 


-0.60 


-2.20 


12.35 + 0.19 


-3.47 + 0.23 



" Fixed slopes of the infrared LF. 



Table A.5. Evolution of the comoving IR energy density and the relative contrib ution of "faint" galax ies (10' Lq < Lir < 10" Lq), LIRGs 
(10' ' Lo < LiR < 10'- Lo) and ULIRGs (Lir > 10'^ Lg). The local reference is taken from 



Sanders et al. 



(20031 



Comoving infrared energy densities 

constrained using our PACS observations are highlighted in bold. Others values rely upon extrapolations based on ^^^e and L^^^e and a faint-end 
slope of the IR LF fixed at its z ~ value. 



Redshift Total IR LD 


"faint" galaxies 


LIRGs 


ULIRGs 


[10*xLoMpc-3] 


[10*'xLoMpc-''] 


[10*xLoMpc-3] 


[lO'-xLo Mpc-3] 



0.00^ 


131!« 


122+^ 


r\ 


50+"-"^ 

^•-'^-0.05 


0.25 


141 + 18 
^^^-18 


117+14 

'-'-'-14 






0.55 


421+85 


220!^'° 


180+40 




0.85 


544+139 


242+^f 




38- 


1.15 


7=9+254 
'-^^-273 




070+143 

■^'^-144 




1.55 


789!2V6 


229!iio 


-350+159 
-'-'^-148 


198!^^ 


2.05 




215+226 


336^354 


538!}^J 
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